THE 


PHYSICAL REVIEW 


A Journal of Experimental and Theoretical Physics Established by E. L. Nichols in 1893 


Vor. 60, No. 5 


SEPTEMBER 1, 1941 


SECOND SERIES 


Scattering of Protons by Magnesium and Aluminum 


T. R. Wi_krns* 
University of Rochester, Rochester, New York 
(July 14, 1941) 


A new type of nuclear scattering camera utilizing photographic emulsions to determine 
number and energy of scattered protons at many sharply defined angles is described and com- 
pared with earlier experimental arrangements. A survey of the effects of various factors in 
producing straggling of the scattered particles is presented. The angular distributions of 
elastically and inelastically scattered proton groups are presented, and the former compared 
with Rutherford scattering. Excited states of magnesium nuclei are found to occur at 1.37 Mev, 


2.80 Mev, and 4.07 Mev. 


INTRODUCTION 


T has been recognized that much valuable 
experimental information concerning the 
fields of force about atomic nuclei of the lighter 
elements can be gained by a direct study of the 
scattering of protons, deuterons and alpha- 
particles by such nuclei.! Yet relatively few 
experimental results have been published in this 
field in recent years,?-!® and the majority of 
these were obtained only with alpha-particles. 


* Because of the untimely death of Professor T. R. 
Wilkins this work has been completed and presented by 
Gerald A. Wrenshall, McMaster University, Hamilton, 
Canada. 
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Published results have appeared in two general 
forms. In one procedure the ratio of observed to 
expected or classical yield of particles scattered 
through a selected scattering angle has been 
measured as a function of the energy of the 
bombarding particles. Another procedure has 
included measurement of the same ratio as a 
function of the angle of scattering for bombarding 
particles of constant incident energy. The 
expected or classical yield is that predicted by 
Rutherford’s theory'® for scattering of particles 
by an inverse square field of force. Where 
the scattering of large numbers of particles 
has been studied, no publications have, until 
recently,®!°"\"4 contained measurements of the 
energy distribution of the scattered particles. 

"4Powell, May, Chadwick, and Pickavance, Nature 145, 
893 (1940). 
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The general lack of intense sources of protons, 
deuterons and alpha-particles until recent years, 
and of suitable means of faithfully registering 
scattered particles for long periods of time, 
has forced investigators in this field to use poorly 
defined scattering angles in order to obtain 
results. This lack of precision has rendered 
theoretical analyses of experimental results 
somewhat inconclusive.'? Both of these de- 
ficiencies have been improved upon in the 
present analysis, in which measurement of the 
yield of scattered particles at a series of sharply 
defined angles is attempted. A cyclotron is used 
to produce a sufficiently intense and homo- 
geneous source of bombarding particles (protons). 
The yield and energy of scattered protons at the 
selected scattering angles are measured from 
the number and length of the discrete tracks 
produced in suitable photographic emulsions by 
the passage of these particles. 


EXPERIMENTAL 


A scattering camera consisting of a flat 
cylindrical container made of brass 0.9 cm thick, 
and of about 30 cm outside diameter and 5.0 
cm depth has been employed. It contains a 


” P. Seligmann, Phys. Rev. 58, 492 (1940). 
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series of sixty-six photographic plateholders set 
fan-wise at five-degree intervals about a central 
scattering foil F (Fig. 1(a)) so that angular 
scattering over the range 15° to 170° can be 
studied. The plateholders are arranged so that 
particles entering the camera at £ and scattered 
by F enter the emulsion surface of each (1’’ 3’) 
photographic plate at a grazing angle of 4°. 
The ensemble of plateholders is milled from a 
single disk of brass. Each photographic plate is 
held firmly in position by a flat bow spring B. 
Between scattering foil and photographic 
plates there is a brass ring R in which are drilled 
the series of 0.09-cm diameter beveled holes 
which, with the scattering foil and the apertures 
through which the incident beam enters, define 
the angles through which the scattered particles 
pass to the photographic plates. The particles 
thus scattered produce flat parallel tracks in the 
photographic emulsion distributed over a very 
restricted area (1.00.1 cm), so that tracks 
produced by extraneous particles are easily 
recognized and eliminated from records. 
Shutters S make it possible to block any of 
the apertures in R during a run if exposures of 
differing amounts are desired for different angles. 
Varying exposures are clearly very necessary, 
especially at smaller scattering angles, if the 
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Fic. 2. Degree of me 


angular definition of 
scattering camera 
shown at two charac- = 
teristic angles. Repre- 
sentative photographic 
plates P are shown. 


distribution of scattered particles varies in ap- 
proximate accordance with Rutherford’s theory. 
Each shutter consists of a brass plate fastened 
to a tapered rotatable brass plug situated in the 
removable side of the camera, and serves to 
cover two adjacent apertures in R. To measure 
exposures, a current integrator at J, attached 
to an ionization chamber registers the beam 
after it has passed through the camera, so that 
suitable relative exposures of the photographic 
plates can be measured. A vacuum is maintained 
throughout the entire camera system by means 
of a pumping system attached at V. 

Stopping foils of aluminum inserted between 
scattering foil and photographic plates reduce 
the ranges of scattered particles by known 
amounts to lengths suitable for registration in 
the plates, and also reduce an effect which is 
roughly equivalent to additional straggling 
which shall be referred to as apparent straggling, 
caused by the inequality in stopping power of 
the component parts of the photographic emul- 
sion. The stopping foils are located in a rec- 
tangular ring trough cut on the inner side of R. 
The scattering foil F is placed at 45° to the beam, 
allowing simultaneous exposure of thirty-three 
plates distributed over diametrically opposed 
forward and backward quadrants of the camera. 
F can be rotated out of the beam, and accurately 
returned to its original position by a tapered 
brass arm A by which it is supported. A extends 
through the back of the camera, and can be 
adjusted during an experimental run. It may be 
replaced during a run by a second similar 
scattering foil arrangement not shown in Fig. 


1(a) by which the remaining two quadrants of 
the camera can be used. 

The angular definition of the camera proper 
is shown in Fig. 2 in terms of the camera’s 
geometry. The individual scattering angles thus 
defined cannot vary from their mean value by 
more than about +1.2°. However, in addition 
to this definition, the beam is further restricted 
as it leaves the cyclotron by two 1.0 mm wide 
slits S (Fig. 1(b)) spaced 16 cm apart, after 
which it passes through an evacuated sylphon 
tube for a distance of 122 cm before reaching 
the camera. If it be assumed that the paths of 
the particles in this section are rectilinear, their 
maximum angular spread on scattering to a 
given photographic plate is reduced to about 
+0.5°. From these considerations it is felt that 
the maximum angular variation does not 
exceed +1°. 

Throughout the present series of experiments, 
a single lot of Eastman Kodak ‘“‘fine-grain 
alpha-particle” photographic plates have been 
used. The relative stopping power of these plates 
was obtained primarily by measuring the lengths 
of alpha-particle tracks of known energy. The 
emulsion sensitization, processing, and micro- 
scope techniques used have been described 
elsewhere."* However, certain refinements have 
been made. To permit systematic counting of 
tracks and to have a check on the magnification 
of the microscope, each plate was marked 
photographically with a series of straight: lines 
perpendicular to the direction of the emulsion 
tracks, and spaced 0.25 mm apart. These lines 


% T. R. Wilkins, J. App. Phys. 11, 35 (1940). 
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were prevented from penetrating deeply into 
the body of the emulsion by making the neces- 
sary exposure with ultraviolet light. This permits 
emulsion tracks located deep in the emulsion to 
be resolved more readily than if more penetrating 
illumination were used. The necessity for using 
depth focusing in emulsion track work is really 
an advantage in that it usually permits an 
observer to distinguish between the ever-present 
“fog” particles in the emulsion and the true ends 
of particle tracks. However, it makes impractical 
the publication of a clear photograph of emulsion 
tracks in this experiment. 

Because of the considerable eyestrain and 
fatigue caused by long periods of close examina- 
tion with a microscope, the microscopic image 
of the tracks was projected onto a_ vertical 
ground glass screen, with ground side toward 
the microscope. In contact with this side was 
mounted a frame carrying a series of twelve fine 
taut steel wires, spaced 0.8 em apart, which 
crossed the track images at right angles and 
permitted the accurate measurement of emulsion 
track lengths. A large plano-convex lens mounted 
on the observer's side of the ground glass screen 
converges the light so that all parts of a track 
image can be viewed simultaneously. This 
arrangement has been found to be much less 
tiring to operate than the usual binocular or 
monocular microscope with eyepiece _ scale. 
Depth focus and the mechanical stage of the 
microscope are operated by the left hand, leaving 
the right hand free to record measurements. 

The intense illumination required for dark- 
field projection was obtained from an automatic 
feed arc lamp. Best results were obtained by 
using small diameter cored carbons to prevent 
wandering of the are spot, and consequent 
fluctuations in the brightness of the projected 
images. The light beam was heat filtered by 
passage through a water solution of copper 
chloride. 


ACCURACY OF OBSERVATIONS 
(1) Factors affecting width of distribution peaks 


It is of some interest to investigate the relative 
importance of the factors which contribute to 
the breadth of the distribution peaks obtained 
by plotting the frequency against length of 
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emulsion tracks, the peaks representing proton 
groups. The factors will include the normal 
straggling of the protons in the homogeneous 
stopping and scattering foils, apparent straggling, 
and error in measurement of track lengths due 
to least count and grain spacing along tracks in 
the photographic emulsion. Superimposed on 
these will be the effects of inhomogeneity of the 
proton beam as it arrives at the scattering 
camera, and, in the case of magnesium, the 
broadening of the peak produced by elastic 
scattering due to the inherent presence of the 
three important isotopes. In Table I are pre- 
sented measurements made at four representa- 
tive angles on the widths of the distribution 
peaks for 6.9-Mev protons scattered by mag- 
nesium using the scattering camera. All of the 
photographic plates were simultaneously ex- 
posed, and were analyzed by the same observer 
using the same microscope. 

Examination of the scattering arrangements 
of the camera as shown in Fig. 2 shows that the 
mean straggling of a homogeneous beam of 
bombarding particles, scattered back through 
the same surface of foil through which they 
entered, will, to a fair approximation, be greater 
than for particles scattered through the opposite 
side by the air equivalent of the scattering foil. 
This statement refers specifically to the arrange- 
ment used and shown in Fig. 1, namely where 
photographic plates are exposed only in the 
diametrically opposed quadrants of the camera 
indicated by the normals to the scattering foil. 
For a homogeneous proton beam there would 
thus be a relatively large increase in the strag- 
gling of scattered particles for scattering angles 
greater than 90°, roughly proportional to 
scattering foil thickness. But when the other 
factors previously mentioned produce distribu- 
tion curves of considerable average width, this 
effect is reduced in importance since it is not 
additive. For example, in this experiment the 


TABLE I. Mean peak half-widths, in cm of equivalent air, of 
distribution curves produced by protons scattered 
elastically (E), and inelastically 
(11) by magnesium. 


125° 145° Average 
2.26 2.57 2.38 
2.07 2.15 2.00 


Scattering angle 60° = 80° 
E peak half-width 2.32 2.35 
T, peak half-width 1.69 2.04 
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average increase in straggling for scattering 
angles greater than 90° would have been +0.62 
cm air for a homogeneous beam of bombarding 
protons, while actually for the representative 
results of Table I it produced a mean peak width 
increase of only +(2.26—2.10)=+0.16 cm air 
straggling. Hence it would appear that this 
factor in broadening the large-angle scattering 
peaks is relatively ineffective in this experiment 
due to the presence of more important factors. 

Because of the dependence of recoil energy of 
a scattering nucleus on-its mass, a beam of 
homogeneous protons scattered by a mixture of 
the Mg isotopes would appear in the elastic 
scatter peak as three components slightly dis- 
placed relative to each other. A simple calcula- 
tion shows that the displacement of the three 
isotope peaks would be slight in the present case, 
being of the order of 0.02 Mev. This would be 
scarcely detectable keeping in mind that Mg* 
constitutes 80 percent of normal magnesium. 
It would at most produce a slight distortion of 
the elastic scatter peak on the high energy side, 
an effect which is not observed here. 

An indication of the apparent straggling 
effect which the inhomogeneous photographic 
emulsion produces can be found by comparing 
the average widths of the two main scatter peaks 
in Mg. The elastic scatter peak should be the 
broader of the two due to the greater path of 
its particles through the photographic emulsion. 
This is seen to be true for the representative 
results in Table I. There is an average increase 
in the width of the elastic scatter peak over that 
of the inelastic scatter peak of +(2.38—2.00) 
= +0.38 cm air. The corresponding increase in 
mean peak width for passage of the proton beam 
through a thickness of air equivalent in stopping 
power to the emulsion thickness between the 
peaks considered is +0.19 cm air. Thus the 
photographic emulsion inhomogeneity roughly 
doubles the straggling in our work." 

Uncertainties in measuring track lengths due 
to the finite spacing of silver grains defining the 
tracks, and by least count error in estimating 
track lengths, are of the same order of magnitude. 
The average grain spacing for protons was found 


% This measurement should be independent of all such 
factors as initial inhomogeneity of scattered beam, and 
characteristic errors in measuring track lengths. 
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NUMBER OF TRACKS 


| 
4 8 12 
LENGTH OF TRACKS IN EMULSION 


Fic. 3. Range distribution curves for protons scattered 
by magnesium through various angles ¢. The curves have 
been shifted laterally to align the two main peaks, a slight 
shift being necessary since the recoil energy imparted to the 
scattering nuclei varies with the angle of scatter. One 
scale division = 19.6u. 


to be about 0.1 scale division in track length. 
Track lengths were measured to the nearest 
tenth of a scale division. It is estimated that, 
combined, these effects produce an average 
uncertainty in track range measurements of 
+0.6 cm of equivalent air. 

Since protons were allowed to enter the 
photographic emulsion at a mean grazing angle 
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Fic. 4. Relative yield of elastically and inelastically 
scattered protons as a function of angle for magnesium 
(run No. ITI). 


of 4°, it is important to know surface and depth 
conditions in the emulsion. Interference fringes 
were produced between the surfaces of samples 
of undeveloped emulsion plates and optically 
flat glass, using sodium yellow light. It was found 
that the emulsion surfaces were surprisingly flat 
except in the immediate vicinity of the natural 
edges of the plates. Edges along which glass cuts 
had been made showed no abnormal effects. 
From plate to plate the emulsion thickness 
varied between thickness of 30u to 554. At the 
4° entrance angle this permits corresponding 
tracks of maximum length between 400 and 750u. 
The finite thickness limits the energy of the 
protons which can be successfully studied by 
this method, although considerably thicker 
emulsion plates are being developed for this 
purpose. 

While it is difficult to estimate the homo- 
geneity of the incident cyclotron beam from 
available experimental data, a minimum for the 
mean energy variation can be tentatively con- 
sidered as +0.1 Mev, which is twice the peak 
voltage across the dees. This is undoubtedly a 
smaller variation than actually was obtained in 
the experiments, since variations in the field 
strength of the cyclotron magnet over the time 
required for a run are not allowed for. 


(2) Precautions relative to scattering angle 


To obtain accurately known and reproducible 
scattering angles care was taken to insure that 
the center of the ellipse defined by the inter- 


section of the beam of incident protons and the 
scattering foil was located at the geometrical 
center of the camera. To prevent reduction in 
the accuracy of scattering angles due to the 
effect of the cyclotron’s magnetic field on the 
scattered protons, the camera was located where 
the field was found to be effectively zero. As an 
added precaution, the plane of scattering of the 
camera was vertical, so that in any case any 
vertical magnetic field could not alter the 
direction of scattered particles in the plane of 
scattering. 


ADVANTAGES AND DISADVANTAGES OF THE 
EMULSION TRACK TECHNIQUE 


An excellent general survey of the basic 
advantages and limitations in the use of tracks 
of nuclear particles in photographic emulsions 
has been recently presented by Shapiro.?° 
Powell and Fertel*! have demonstrated some of 
the marked advantages of the technique in the 
measurement of neutron energies. The undesir- 
able effects of apparent straggling in emulsion 
track studies have been discussed by Taylor.” 

Outstanding advantages of the emulsion track 
technique in nuclear scattering experiments are: 
(1) Simultaneous exposures can be made at a 
large number of sharply defined angles in a 
minimum of time. (2) Both yield and energy 
distribution of scattered particles are obtained 
at each scattering angle with the simplest of 
apparatus. (3) Since the tracks can be seen, 
random “background” tracks such as_ those 
caused near a cyclotron by neutrons, can be 
recognized as such. (4) Scattered protons of 
energy as low as 0.8 Mev can be readily detected 
and their energy estimated. Nothing in the 
nature of an entrance window is required for 
emulsions, as in the case of most impulse 
counters and cloud chambers. 

Because of the comparatively slow rate at 
which track lengths can be measured,” several 
observers have taken part in the track analysis, 
using three independent microscopes set to the 
same magnification. In order to make a check 


20 M. Shapiro, Rev. Mod. Phys. 13, 58 (1941). 

2t Powell and Fertel, Nature 144, 115 (1939). 

2H. J. Taylor, Proc. Roy. Soc. A150, 382 (1935). 

*3 About 200 tracks per hour can be accurately measured 
for a limited time. 
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on the personal error involved in such mixed 
analyses, two observers, working independently 
on different analyzers, in turn analyzed the same 
plate (scattering angle 145°). A comparison of 
these two analyses is included in Fig. 3 in which 
it is seen that no important variation in curve 
shape present in the one is absent from the other. 
The slight difference in total track number is 
believed due to the fact that one observer 
counted only those tracks appearing in the 
portion of the plate seen by the scattering foil, 
while the other examined the entire plate surface 
for additional stray tracks. 


RESULTS AND DISCUSSION 


(1) Yield of scattered protons as a function of 
scattering angle 


As a preliminary test of the scattering camera, 
6.6-Mev protons were scattered by 0.2-mil 
platinum foil (run No. I), and the relative 
numbers (yield) scattered through five different 
angles were measured. The variation in yield 
with scattering angle agreed with that predicted 
for Rutherford scattering, as was expected. 
Subsequently, the scattering of 6.9-Mev protons 
by 0.5-mil foils of aluminum and magnesium 
has been studied and preliminary results pre- 
sented® !%'* announcing the detection of inelastic 
scattering for both elements, and the variation in 
the relative amounts of inelastic and elastic scat- 
tering with angle for the latter. 

The apparent straggling effect produced in 
the photographic emulsion was reduced in the 
aluminum run (No. II) and in one of the mag- 
nesium runs (No. III) by insertion of a suitable 
thickness of aluminum stopping foils between 
scattering foil and photographic plates. This had 
the effect of increasing the resolution between 
the higher energy groups of scattered protons 
while cutting off all scattered protons of less 
than 2.5-Mev energy from the photographic 
plates. It also shortened the resulting emulsion 
tracks to such a degree that they could be 
conveniently examined as units under a micro- 
scope using a 10 X eyepiece and a 4-mm objective. 

Track numbers and lengths have been meas- 
ured for seven selected plates of run III, and 
the series of distribution curves so obtained is 
shown in Fig. 3. The peaks appearing at 12.2 scale 
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divisions are believed to be produced by elasti- 
cally scattered protons. This conclusion is based 
on a comparison of the energy of the group of 
scattered protons of highest energy, and the 
energy of the incident ones. By use of range 
energy graphs** and the measured lengths of 
emulsion tracks, the former were found to have 
a mean energy of 6.87 Mev minus nuclear recoil 
energy, while from air range measurements the 
latter were found to have 6.9-Mev energy. The 
peak at 6.9 scale divisions is seen to appear in 
all plates of run No. III which have been 
analyzed so far, and the marked variation with 
angle of its intensity NV, relative to the intensity 
Ne of the peak at 9.1 scale divisions is shown 
in Fig. 4. The ratio N;/Ng appears to pass 
through a maximum value for 90° scattering. 
The angular distribution of protons scattered 
by aluminum (run No. II) differs greatly from 
that predicted for Rutherford scattering. (Fig. 5, 
curves B and A, respectively.) While curve B 
includes the total number of protons scattered 
by aluminum having energies between 3.0 and 
6.9 Mev, it actually represents the yield for 
elastic scattering alone except at large scattering 
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Fic. 5. Yield of elastically scattered protons from alumi- 
num (curve B) and magnesium (curve C). Curve A repre- 
sents the distribution predicted for Rutherford scattering, 
and is normalized to curve B. 


* M.S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 245 (1937). 
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Fic. 6. Range distribution of protons scattered by 
magnesium for a scattering angle of 95° (run No. IV). 
One scale division = 39.2.. 


angles, since the one observed group of inelasti- 
cally scattered protons is small relative to the 
elastically scattered group, and appears in this 
run only at scattering angles greater than 145°. 
The curve for magnesium (Fig. 5, curve C) is for 
elastically scattered protons alone. These curves 
demonstrate clearly that, exclusive of nuclear 
collisions in which the scattered particle gives 
excitation energy to the scattering nucleus, 
important deviations from Coulomb forces occur 
in the immediate neighborhood of the nucleus. 

There is some evidence in Fig. 3 of unresolved 
proton groups of comparatively low intensity 
at about 9.5 and 11.0 scale divisions. An irregu- 
larity corresponding to the one at 11.0 scale 
divisions has been observed in a_ different 
experimental study now in progress*® of proton 
scattering by magnesium. 

In order to observe the distribution of emulsion 
track lengths over the largest possible energy 
range, another run (No. IV) was made in which 
no stopping foils were used between the mag- 
nesium scattering foil and the photographic 
plates. This allowed representative scattered 
protons of all available energies to reach the 
photographic emulsions. In Fig. 6 the distribu- 
tion in emulsion track length of protons scattered 
through an angle of 95° and measured at one-half 
of the previous magnification, is shown. The 
two large peaks previously examined in run 


2 R. H. Dicke and J. Marshall, Phys. Rev. 59, 917 
(1941). 


No. III are at 5.37 and 8.35 scale. divisions, 
while two further peaks are observed at 2.9 and 
1.3 scale divisions. The presence of both of these 
peaks has been confirmed from the distribution 
curves of two other plates of the same run. 
Collection of information concerning their in- 
tensity variation with scattering angle is at 
present in progress. 


(2) Excitation energies of magnesium nuclei 


An explanation of the series of peaks obtained 
in Figs. 3 and 6 is that an effect analogous to 
that found in the Franck-Hertz experiment is 
occurring when the bombarding protons collide 
with atomic nuclei. While some of the protons 
undergo elastic nuclear collisions, others lose a 
discreet amount of energy to the nucleus. This 
energy may be considered as exciting the nucleus 
to a higher quantized energy state, from which 
it will subsequently revert to its normal state 
with the emission of gamma-radiation. Hence a 
measure of the decrease in energy suffered by 
the groups of inelastically scattered protons 
should represent directly the excitation energies 
of the struck nuclei; and these excitation 
energies and some of their differences should 
represent possible gamma-ray energies which 
the nuclei may emit. 

The energies of the proton groups were 
calculated by determining the equivalent proton 
air ranges from their measured track lengths, 
and the known stopping powers of the photo- 
graphic emulsion, stopping, and scattering foils. 
It is worthy of note, however, that it is incorrect 
to take the excitation energies as the direct 
differences between the energy of the elastically 
scattered group and the energies of the inelasti- 
cally scattered ones as has been done with earlier 
results.'°'’ A correction must be made since the 
recoil energy loss to the nucleus in the inelastic 
collisions differs from nuclear recoil loss in elastic 
collisions.** Referring to Fig. 6 the proton 
energies corresponding to peaks EF, J;, J2, and J; 
are 5.94 Mev, 4.62 Mev, 3.25 Mev, 2.04 Mev, 
respectively, and the excitation energies corre- 
sponding to peaks J;, J2, J3 are 1.37 Mev, 2.80 
Mev, 4.07 Mev, respectively, in satisfactory 
agreement with values of Dicke and Marshall.* 


% ] wish to thank Mr. Dicke for developing the necessary 
correction equations which were used here. 
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It is only possible with the information 
available to offer a partial indication as to which 
of the proton groups J;, Zs, J; is associated with 
which of the nuclei Mg™, Mg*, It is 
reasonably certain, however, that the group J; 
is produced by Mg*, since at a scattering angle 
of 90° it is nearly 1.7 times as large as the E 
group and could scarcely be produced by the 
other two rarer isotopes. This observation is 
supported by Henderson's” measurement of the 
gamma-ray energy of radio-sodium in the 
proposed reaction, 


Na*>Mg*+8+7, 


in which he finds, by the absorption coefficient 
method, that the gamma-ray energy has an 
approximate value of 1.3 Mev which compares 
with our value of 1.37 Mev. However, from a 


7 M.C. Henderson, Phys. Rev. 48, 855 (1935). 


number of published results, Livingston and 
Bethe* list gamma-rays for this reaction of 
energies 0.95 Mev, 1.93 Mev, and 3.08 Mev, 
none of which is closely comparable with the 
excitation energies or their differences as reported 
here. 

From published studies of the proton groups 
from the reaction 


Na®*+ 


Livingston and Bethe* have assigned excitation 
energies to Mg*® of 2.2 Mev, 4.0 Mev and 5.0 
Mev. Of these, the first mentioned is not ob- 
served in our results, while the last value quoted 
is not in the range of our experiment. However, 
the group J; of measured excitation energy 4.07 
Mev might perhaps correspond to the remaining 
proton group, thus associating this group with 
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The forced separation of variables usually employed in 
the quantum-mechanical treatment of molecular problems 
introduces certain small errors into the wave function. If 
the exact Hamiltonian is used, energy values can be 
computed very accurately because first-order errors in a 
wave function give rise to second-order errors in the energy. 
The energy of H:* is ordinarily computed by using a 
separable, approximate Hamiltonian instead of an exact 
one. From a consideration of the terms which must be 
added to the approximate Hamiltonian to make it exact, 
Van Vleck has derived the correction terms needed to 
reduce the first-order error in the computed energy to a 


1. THE PROBLEM 


HE usual quantum-mechanical treatment 
of molecular problems is based upon a 
forced separation of variables.' This introduces 
certain small errors into the wave function, 
which, as such, is still as accurate as necessary. 


Cf. A. S. Coolidge and H. M. James, J. Chem. Phys. 6, 
730 (1938). 


second-order error. The correction term is a function of R, 
the internuclear distance, and its calculation requires a 
knowledge of the wave function. In this paper the ground 
state wave function of H2* is accurately determined over 
values of R from 1.20 to 2.75 atomic units and a table of the 
wave function coefficients is given, along with the corre- 
sponding energy values. Then, correction terms are 
calculated for a set of values of R. Including the proper 
correction term, the total negative energy of H.* for the 
equilibrium internuclear distance is found to be 1.20472 
+0.00001 Ey = 132,132+10 cm. This result is compared 
with the consequences of certain experimental data. 


Since energy values are subject to much closer 
check by experimental data, it is desirable to 
determine energies as accurately as possible. 
They can be computed very accurately in some 
cases because first-order errors in a wave func- 
tion give rise to second-order errors in the energy 
if the exact Hamiltonian is used in calculating 
the energy. The energy of H2* is computed in 
the ordinary way by using, not the exact Hamil- 
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tonian, but the separable, approximate one. 
From a consideration of the terms which must 
be added to the approximate Hamiltonian to 
make it exact, Van Vleck? has derived the cor- 
rection terms needed to reduce the first-order 
error in the computed energy due to total neglect 
of the motion of the nuclei to a second-order 
error due to carrying out the computation with 
approximate wave functions. The correction 
term, denoted by Prn+QnatSnn, is a function 
of R, the internuclear distance. The potential 
energy at each value of R is obtained by adding 
this correction to the energy value of the corre- 
sponding fixed-nuclei solution. For infinite sepa- 
ration this correction becomes that required to 
convert the fixed- into the free-nuclei problem, 
hence it may be termed ‘“‘the correction for 
nuclear motion.” In this paper the ground state 
wave function of H.* is accurately determined 
over values of R from 1.20 to 2.75 atomic units 
and the corresponding corrections for nuclear 
motion evaluated. A corrected value of the 
energy of H.*+ for equilibrium internuclear dis- 
tance is obtained and compared with certain 
experimental results. 


2. CALCULATION OF THE WAVE FUNCTION 


The calculation of the wave functions of H.* 
has been discussed in numerous papers in recent 
years;*-"' but, due to the difficulties of calcula- 
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tion described below, the various authors have 
calculated only the energy values and not the 
wave function itself. Hence the first step in 
calculating the correction for nuclear motion is 
to make an accurate determination of the ground 
state wave function. 

If elliptical coordinates, £, 7, y, are inserted 
in the Schrédinger equation, the latter becomes 
separable and the wave function for the ground 
state may be written: 


Cc 
(27) 


where C is the normalization constant. Y() and 
X(£) are defined by the differential equations: 


d dY 
+n!) Y=0, (2) 
dy 


d dX | 
—(1— + 2REF =0, (3) 
dé dé 


in which \?= — E’R?/2, E’ being the electronic 
part of the energy. 
The solution of (2) is given by 


Y() =1+42(A)P2(m) +a4(A) P4(m) 
+a6(A)Pe(n) +--+ (4) 


where, in order that the series may converge, the a; must satisfy the relation: 


oA? +a 
(21—3)(21—1) 


(41)? 
(21-41)(21+3)  (214+1)(2/—1) 


=0. 
(214+-3)(2/+5) 


So that the above equation may be satisfied by values of a; other than zero, uv’ and \? must satisfy a 
functional relation, given in continued fraction form by Wilson‘ and solved in series form (accurate 


to the sixth decimal place) by Sandeman:® 


92,” 


2894282\"4 


(5) 


2 J. H. Van Vleck, J. Chem. Phys. 4, 327 (1936). 

3Q. Burrau, Kgl. Danske, Vid. Selsk. 7, No. 14 (1927). 

* A. H. Wilson, Proc. Roy. Soc. A118, 617, 637 (1928). 

5 P. M. Morse and E. C..G. Stueckelberg, Phys. Rev. 33, 
932 (1929). 

°E. Teller, Zeits. f. Physik 61, 458 (1930). 


3 33-5 35-5-7 39-5?-7-11 


7E. A. Hylleraas, Zeits. f. Physik 71, 739 (1931). 

5G. Jaffe, Zeits. f. Physik 87, 535 (1934). 

® I. Sandeman, Proc. Roy. Soc. Edinburgh 55, 72 (1935). 

10 W. G. Baber and H. R. Hasse, Proc. Camb. Phil. Soc. 
31, 564 (1935). 

1S. K. Chakravarty, Phil. Mag. 28, 423 (1939). 
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TABLE I. Wave function coefficients and energy values. 


R as a4 a6 by be bs E* E* 
1.0 1.214301 0.114518 0.001978 0.0000137 0.008496 0.0003706 0.0000411 3.612652 —2.712737 —1.065698 
11 1.366077 -139388 002919 -.0000243 010137 .0003908 -0000406 3.660361 —2.593552 —1.129506 
1.2 1.523336 -166929 004159 -0000410 011831 .0004032 .0000392 3.723734 —2.482167 —1.169259 
1.3 1.685909 .197207 .005770 .0000663 013559 -0004086 .0000373 3.802924 —2.378365 —1.192061 
1.4 1.853595 .230283 007815 0001035 015301 .0004077 0000348 3.898359 —2.281847 —1.202863 
1.5 2.026162 .266212 010364 0001565 .017039 0004014 .00003 23 4.010611 —2.192272 —1.205184 
1.6 2.203349 .305038 013487 .0002306 018757 .0003906 0000294 4.140736 —2.109276 —1.201567 
1.7 2.384868 346787 017253 .0003317 020438 0003762 0000267 4.292105 — 2.032493 —1,.193872 
1.8 2.570407 391465 021719 0004674 022068 -0003590 -0000239 4.459978 —1.961558 —1.183471 
1.9 2.759633 439049 .026927 -0006465 023634 .0003399 0000214 4.653078 —1.896115 —1.171381 
1.483559 1.997466** 0.260107 0.009908 0.0001465 0.016754 0.0004028 0.0000328 3.991018 —2.206537 —1.205268 


* FE’ and E are given in terms of Ey =109,679 cm =13.530 ev. 
** Equilibrium internuclear distance. 


The solution of (3) is written 


{1 +0400) —) (6) 
E+1 
where p=1—(R X). In order that the series may converge, },, must satisfy the relation 
= — Unda Dn, ly by =0, (7) 
where 
+’) 
— =—— ))*, (m+1)?*. 
(m+1)* 


So that (7) may be satisfied by values of },, other than zero, yn’, A, and p (or R) must be related in the 
manner indicated by the continued fraction: 
Vo Veo 
—u=— -- — (8) 
Ur 


In theory it is possible to eliminate yp’ between (5) and (8) to obtain a relation between \ and 
R;'°" in practice it is extremely difficult to determine an algebraic relation between \ and R suffi- 
ciently accurate to yield reliable values for the wave function coefficients. In the present work this 
calculation has been done numerically to insure maximum accuracy; for an arbitrary value of \, »’ 
was found from (5), and then (8) was solved for the corresponding unique value of p (and hence R). 
From these values (all accurate to the sixth decimal place), corresponding values of a), bn, C, and the 
energies E’ and E[ = E’+(e?/R)] were determined with like accuracy. Table I gives these results 
over a range of values of \ (or R); the last row gives results for the equilibrium value of R. 

The value of the energy for the equilibrium value of R is exactly equal to that found by Hylleraas.’ 
Sandeman? calculated a somewhat smaller absolute value. His result probably is slightly in error 
because the empirical formula relating \ and R, from which he obtained the energy minimum by the 
appropriate differentiation, was obtained from pairs of values of \ and R which fail to give agreement 
in Eq. (8) by about two percent. 

The pairs of values of \ and R tabulated in Table I give rise to the empirical relation: 


R=0.0099234 +.0.9302959d +.0.2461946d? + 0.0449322* — 0.0176366A4 +0.0006485\° 
— 0.000097 (9) 


3. CALCULATION OF Pyn+QnatSnn AND DiscUSSION OF RESULTS 


Van Vleck’s® expression for S,.,, which represents the energy difference due to the fact that the 
center of gravity of the molecule does not exactly coincide with the center of gravity of the nuclei, 
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(10) 


where 1/p=mass of proton, Ey, =ionization potential of hydrogen. 
The corresponding expression for Pxn+Qnn, Which may be qualitatively described as representing 
coupling between vibrational and electronic motion, is 
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The expression defined in Eq. (10) has been 
evaluated in a straightforward manner for a 
series of values of \. The evaluation of expression 
(11) involves the calculation of the derivatives of 
the coefficients a; and b,, and of \ and p with 
respect to R." The appropriate derivatives have 
been calculated for different values of \ by use 
of the Gregory-Newton interpolation formula 
for derivatives of tabulated functions. Otherwise 
the evaluation of (11) involves only lengthy, 
but straightforward, integrations. Figure 1 shows 
the variation of PantQnn, and 
+S,, with internuclear distance. The total 
correction for nuclear motion at the equilibrium 
distance is 0.00055 +0.00001 E,,(60 cm~! = 0.0075 
ev). Van Vleck! has estimated this quantity as 
114 cm“'=0.014 ev. 

By combining the correction for nuclear mo- 
tion with the value of the energy at equilibrium 


” The normalization constant C is also a function of R, 
but calculation shows that terms involving 8C/@R are 
negligible compared to the others; hence such terms are 
omitted in Eq. (11). 


given in Table I, the total 
negative energy of H:*+ becomes, for this value 
of R, 


= 1.20472 +0.00001 Ey, 
= 132,132+10 cm™ 
= 16.300+0.001 ev. 


(12) 


This value now can be compared with the results 
of.experiment. 

The total energy of He equals twice the 
ionization energy of the H atom plus the dissocia- 
tion energy of He. Beutler and Junger™ have 
obtained a direct experimental value of the 
dissociation energy: 


= 36,116+6 cm". 
With this value, the energy of Hy. becomes 


W (He) = 2(109,679) + 36,116 
= 255,474+6 cm". 


(13) 


18H. Beutler and H. O. Junger, Zeits. f. Physik 101, 304 
(1936). 
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Fic. 1. Variation of Sn» and 
+Qnn With internuclear distance. 


The ionization energy of Hz can now be calcu- 
lated, for the total energy of He equals the 
ionization energy of He plus the energy of H.2* 
minus the zero-point vibrational energy of He*. 
By extrapolation from a singlet series of He, 
Richardson’ has obtained for the zero-point 
energy the value 0.1399+0.0007 ev; Sandeman’s 
theoretical work gives 0.1390 ev. The computa- 
tions of this paper lead to the value 


0.1389 ev=1126 
(error not exceeding 10 cm='). (14) 


4 QO. W. Richardson, Proc. Roy. Soc. A152, 503 (1935); 
Molecular Hydrogen and Its Spectrum (Yale University 
Press, New York, 1934), p. 158. 


This value of the zero-point energy and the 
value of the energy of H,* given in (12) give for 
the ionization energy of He: 


I(He) = 255,474 — 132,132+1126 (15) 
= 124,468+ 20 


Beutler and Junger, noting that the ionization 
energy of H» cannot be measured with the same 
accuracy as the dissociation energy, have ob- 
tained the experimental value 


1(H2) = 124,427+10 


The small discrepancy may indicate that the 
probable error has been underestimated in one 
or both cases. 

The dissociation energy of H+ equals the 
energy of H2* minus the zero-point vibrational 
energy of H.* minus the ionization energy of the 
H atom. Using the theoretical values (12) and 
(14), one obtains: 


= 132,132 — 1126— 109,679 (16) 
= 21,327+20 


This result cannot be directly compared with 
experiment because direct observations on the 
H.* spectrum are completely lacking. The 
experimental value, 21,366+15 cm™, given by 
Beutler and Junger is a combination of the 
ionization and dissociation energies of He; hence 
the small discrepancy between calculated and 
experimental values of J(H2) appears again here. 

The author wants to thank Dr. Hubert M. 
James for suggesting this problem and for many 
discussions concerning its solution and Dr. K. 
Lark-Horovitz for his encouragement during 
the progress of the work. 


| | 
| | 
: | 
| | 
\ | | | 
\ | | | 
| 
Pan + Qnn | 
R 
1) 
al 
1e 
2) 
1¢ 
ve 
_| 


SEPTEMBER 1, 1941 


PHYSICAL REVIEW 


VOLUME 60 


The Multiple Scattering of Fast Electrons 


N. L. Ovreson,* K. T. Cuaot ann H. R. 
The Harrison M. Randall Laboratory, University of Michigan, Ann Arbor, Michigan 


(Received July 11, 1941) 


An extensive cloud-chamber investigation of the scattering of fast electrons (3-11 Mev) by 
foils of various elements has been completed. The foils were of such thickness that the resultant 
average deflection was due chiefly to many deflections through small angles. This fact makes 
possible a comparison with two recent theories of multiple scattering. The results indicate that, 
for elements at either end of the atomic number scale, the scattering is less than predicted by 
theory. For elements of intermediate atomic numbers, there is good agreement with theory. An 
attempt was made to detect the effect of crystal size on scattering. The results were negative 
within our limits of error, and, if such an effect exists, it must be small. 


INTRODUCTION 


NE important aspect of the problem of 

scattering of electrons is the scattering 
through extremely small angles. This usually 
means scattering in which the electrons pass 
through that part of the atomic field which is 
partially screened by the orbital electrons. So 
far, we have no method of studying this kind of 
scattering as individual events, because the 
deflections are far too small. Instead, we must 
attack the problem by observing the result of 
multiple deflections. The usual method is to 
allow electrons of known energy to pass through 
a lamina of solid material and to observe the 
difference in direction of incidence and emergence 
for each electron; the angle measured being in 
each case the result of thousands of individual 
deflections. The theory of multiple scattering 
begins with an assumption as to the elementary 
scattering law (which is complicated by the 
screening) and from that point on the problem 
is one of statistically summing the many de- 
flections. Naturally, it is only with the aid of 
an accurate treatment of this statistical problem 
that a comparison of the experimental results 
with theory can be considered to be a good test 
of the laws of interaction between atoms and 
rapidly moving electrons at large distances from 
the nucleus. Recently, two theories’? of multiple 


* Now at the United States Coast Guard Academy, New 
London, Connecticut. 

+ Present Address: Institute of Radium, National 
Academy of Peiping, Shanghai, China. 

‘E. J. Williams, Proc. Roy. Soc. 169, 531 (1939); Phys. 
Rev. 58, 292 (1940). 

2S. Goudsmit and J. L. Saunderson, Phys. Rev. 57, 24, 
552 (1940); Phys. Rev. 58, 36 (1940). 


scattering have been published in which the 
statistics has been sufficiently developed to 
allow a significant comparison of theory with 
experiment. These will be compared with our 
experimental results. 


EXPERIMENTAL PROCEDURE 


Apparatus 


The source of electrons was Li’, obtained by 
the reaction 


Li?+ 


The Li’ disintegrates with a half-life of 0.8 
second according to the reaction 


Li’—Be*+e- 


giving a beta-ray spectrum whose upper limit 
is 12 Mev.? 

The deuterons were accelerated by means of 
the high voltage vacuum tube apparatus, which 
has been described elsewhere.‘ Because of the 
short half-life it was necessary to prepare the Li*® 
just before each expansion of the cloud chamber.*® 
A target of lithium metal was bombarded with 
0.8-Mev deuterons for 2 seconds, then about 
} second was allowed for the cloud chamber to 
be cleared of ions before the expansion occurred. 
The high voltage tube and the cloud chamber 
were controlled by the same motor-driven 
contact system, so that the cycle repeated 
automatically. The electrons emerged from the 


3D. S. Bayley and H. R. Crane, Phys. Rev. 52, 604 
(1937). 

4H. R. Crane, Phys. Rev. 52, 11 (1937). 

5H. R. Crane, Rev. Sci. Inst. 8, 440 (1937). 


378 


= 


An 


SCATTERING OF FAST ELECTRONS 379 


vacuum tube through a thin aluminum window, 
and were allowed to enter the cloud chamber 
through another aluminum window in its side 
wall, falling perpendicularly upon the scattering 
lamina, as shown in Fig. 1. The electrons were 
roughly selected as to energy by the appropriate 
placing of lead slits between the tube and the 
cloud-chamber window, making use of that part 
of the magnetic field which extends outside the 
chamber. To reduce the possibility of confusion, 
the intensity of the electron source was adjusted 
so that, on the average, only about three tracks 
appeared on each photograph. 

For use later in considering possible errors 
due to geometry, it is necessary to state the 
exact dimension of the cloud chamber. It is used 
in a horizontal position, is 15 cm in diameter, 
and has a total depth of 5 cm when expanded. 
The collimated light beam illuminates a slice of 
the chamber 2.5 cm in depth, located in the 
center with respect to the top plate and floor of 
the chamber. The scattering lamina was attached 
to the top glass plate of the chamber, and 
extended well beyond the illuminated part of 
the chamber, so that there was no possibility 
of photographing tracks of electrons which 
passed under the lamina. The window through 
which the electrons entered was 8 mm high and 
2.5 cm long, and was at a distance of about 30 
cm from the lithium source. This means that the 
incident electron paths were very nearly parallel 
to the plane of the chamber, and that they 
struck the scattering lamina at its center. 
Photographs were taken from a point directly 
above the chamber with a Sept 35 mm camera, 
in the vertical position. 


Criteria for selecting and measuring tracks 


Only those tracks were included in the data 
which satisfied the following criteria: (1) The 
incident track must make an angle not greater 
than 15 degrees with the normal to the scatterer, 
measured in the plane projection. (2) The 
incident track must be clear and sharp over a 
distance of at least 5 cm and to within 5 mm of 
the scatterer. (3) The emergent track must have 
a visible length of not less than 7 cm and extend 
to within 5 mm of the scatterer. The illuminated 
portion of the chamber is approximately 2.5 cm 
deep and the point of scattering is located at the 


DEUTERON BEAM 


ALUMINUM 


INCIDENT TRACK 
™ 


SCATTERING 
FOIL 


SCATTERED TRACK \ 


Fic. 1. Arrangement of apparatus. 


middle of the illuminated region. (See Fig. 1.) 
This means that the vertical component of the 
actual scattering is limited to 2.5/2K7=10 
degrees. (4) The track must not be visibly 
distorted by scattering in the gas of the cloud 
chamber. 


Measurement of tracks 


Because of the large number of scattering 
angles which were to be measured, we found it 
necessary to devise a rapid method of measure- 
ment. A set of circles, in steps of } cm, was 
engraved on a piece of transparent celluloid 
about 6X8 inches in size. This was then cut in 
the center, on a common radius of the circles. 
To measure the incident part of the track, a 
circle was fitted to the projected image track as 
shown in Fig. 2. The radius of curvature was 
recorded, and a line drawn along the edge of the 
celluloid at AA. This line is perpendicular to 
the track at the scattering foil. A similar meas- 
urement of the other half of the track is made 
with the other half of the celluloid, drawing a 
line along edge BB. The angle between the two 
lines gives the angle of scattering; also both 
radii of curvature are obtained. The average of 
the two radii is used, and this automatically 
takes care of the energy loss within the scattering 
material. This loss is appreciable only in the 
carbon scatterer, however. 


Scattering foils 


The scattering foil in each case was made of 
such a thickness that the average angle of 
scattering was roughly 10 degrees. A greater 
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thickness (and greater average angle of scat- 
tering) would have made the errors due to 
geometrical limitations of the cloud chamber 
too great. On the other hand, thinner scatterers 
than those used would have given rise to large 
relative errors in the measurements of the angles. 
All the scattering foils used had, therefore, about 
the same scattering power, or, the number of 
atoms per sq. cm area times Z* was about 
constant. 


THEORY 


Full descriptions of the several theories will 
not be attempted here, but only enough of each 
will be given to enable the reader to understand 
clearly how we have used the theory. 


The Williams theory 


In the multiple scattering theory of E. J- 
Williams, the distribution of the horizontal 
projection of the resultant angle of scattering 
(angle @ in Fig. 1) due to individual deflections 
between ¢;, and @min is taken to be Gaussian. 

¢, is defined as that angle such that, on the 
average, each electron experiences one deflection 
greater than @¢; in passing through the foil. 
gmin is an effective cut-off angle, due to the 
shielding of the nuclear field by the orbital 
electrons, the value of which is obtained by 
using the Fermi-Thomas distribution of electrons 
in the atom. A correction for the effect of 
individual deflections greater than ¢; is applied 
later, and its effect on the distribution curve is 
to broaden the Gaussian part and to add a tail, 
which is usually referred to as the “single 


Fic. 2. Method of 
measuring curva- 
ture of tracks. 


scattering tail.” Since the average angle of 
scattering, for a monochromatic electron beam, 
varies inversely as //p, it is convenient to make 
plots in which the abscissa is //pa, as in Figs. 3 
to 6. The angle of scattering of each electron is 
multiplied by its own //p value. In this way, a 
group of electrons having a wide range in J/p 
may be placed on the same plot, for the purpose 
of comparison with theory. 
Williams’ formula is 


= 0.88A $[(9.41+logA —0.67 logZ)! 
+1.81 ]-10° gauss centimeters degrees, (1) 


where A is a parameter characterizing the 
scatterer given by 


A=0.1510Z?/ M, (2) 


where .V/ is the atomic weight of the scatterer of 
atomic number Z. ¢ is the superficial density of 
the scatterer in grams per square centimeter. 
IIp is the product of the magnetic field intensity 
and the average radius of curvature of the 
electron in the field. 


The Goudsmit and Saunderson theory 


The theory of Goudsmit and Saunderson does 
not divide the individual deflections into regions 
above and below a certain angle, such as ¢), in 
the Williams theory, but, by means of standard 
methods of the calculus of probability, the whole 
range is treated at once. (J/pa), is not given by 
a formula analogous to Williams’ formula, 
because the momentum dependence upon @ is 
not linear and explicit. However, the exact 
values of (J7pa),, for a large number of different 
scatterers and energies have been computed by 
Goudsmit and Saunderson,® and, by making 
interpolations of values given in their table, 
(ITpa), is easily obtained for any case which has 
so far been tried experimentally. The scattering 
distribution is essentially Gaussian provided 4& 
is not large (~10 degrees). The corresponding 
average involving the total angle of scattering @ 
(Fig. 1) is obtained by the same procedure. In 
this case, however, the variation with energy is 
found to be small for the range of energies dealt 
with in this investigation and can be represented 


6S. Goudsmit and J. L. Saunderson, Phys. Rev. 57, 24 
(1940), Table III. 
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Fic. 3. Scattering of electrons by C, 0.82 g/cm*. 


to a good approximation by the formula 


= 1.70A 4(5.60—3 logZ+3 logA)!-10° 
gauss centimeters degrees. (3) 


Scattering by orbital electrons 


To account for the scattering due to collisions 
with the orbital electrons, we have assumed that 
each orbital electron scatters as if it were a 
nucleus of atomic number one. Since there are 
Z electrons per atom, we have simply to substi- 
tute Z*+Z for Z* in the theoretical formulas. 
In all the theoretical values quoted by us, this 
substitution has been made. 


The effect of crystal structure upon multiple 
scattering 


The scattering materials used in most of the 
experiments so far have been of polycrystalline 
structure (metals, for example). Also, an im- 
portant part of the multiple scattering is due to 
individual angles of deflection which are of the 
same order of magnitude as the angle between 
maxima in the crystal diffraction pattern for 
electrons of the same energy. Wheeler’ has 
suggested that, in view of this fact, the scattering 
may not be independent of the crystalline 
structure, and that some of the second-order 
effects, such as extinction of the scattered beam 
by rescattering back into the original direction, 
may have an influence upon the total scattering, 


"J. A. Wheeler, Phys. Rev. 57, 352A (1940). 
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Fic. 4. Scattering of electrons by Cd, 0.111 g/cm?. 


as is well known in the case of x-ray diffraction. 
Some preliminary estimates of the magnitude of 
the effect have been made by Wheeler. It is 
interesting to note two ways in which the 
diffraction of electrons differs from that of 
x-rays: (1) Coherence exists only over the 
distance between inelastic collisions, which is 
about 10-5 cm in carbon at this energy, for 
example. (2) The scattering power of each atom 
is so great, that, before many crystal planes 
have been traversed, all the intensity has 
disappeared from the original beam. The partici- 
pation of such a small number of planes in the 
coherent scattering process means that the 
maxima are broad, and the effect of extinction 
is thereby lessened. The first of these limitations 
was pointed out to the authors by Professor 
Wheeler, and the second by Professor Wigner. 


ERRORS AND CORRECTIONS 


Geometrical limitations in the cloud chamber 


It has already been mentioned that a track 
was included in the data, only if the emergent 
portion had a length of at least 7 cm. Since the 
illuminated section of the chamber has a depth 
of 2.5 cm, such a criterion is equivalent to 
imposing an upper limit of about 10 degrees on 
the vertical component of the scattering angle. 
In the region of multiple scattering, where the 
distribution is approximately Gaussian, such a 
limitation would not affect the shape of the 
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z 
Fic. 5. Scattering 
of electrons by Cu 
and Fe. Upper 
curve Cu, 0.174 
g/cm*. Lower curve 
Fe, 0.173 g/cm?. 
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distribution in a, the projected angle of scat- 
tering, whereas, in the region of the single 
scattering tail, this restriction has the effect of 
reducing the single scattering tail and, hence, 
the contribution of single scattering to the total 
scattering. Williams takes account of this 
limitation and gives an approximate formula for 
calculating the reduction in the value of (J7pa),,.8 
Goudsmit and Saunderson correct for this effect 
by a method which amounts to interpolating 
between the theoretical values for two extremes: 
an infinitely thin chamber, and an_ infinitely 
thick chamber.’ Both of these cases can be 
evaluated theoretically, and the difference be- 
tween them is not large. In all the theoretical 


® The formula referred to is Eq. (13) in Williams 
second paper for the case ¥ <¢2, where y is the upper limit 
imposed by the cloud-chamber geometry on the vertical 
angle (~10 degrees) and ¢2, the angle corresponding to the 
intersection of the Gaussian distribution with the single 
scattering curve. (See Fig. 1 of Williams’ first paper) 
¢2>25° for the scattering foils used in this investigation. 

The actual (Hpa),, should lie between and 
(Hpa),, obtained from interpolation of their tables, and 
should be closer to (Hp@)a,. This follows from the fact that, 
if the upper limit on the vertical angle is small, that is, 
only a very thin horizontal slice of the cloud chamber being 
considered, the actual (Hpa)y should be given by (Hp@)ay 
interpolated. On the other hand, if all vertical angles are 
possible, then the actual value is given by (Hpa),y inter- 
polated. As the upper limit on the vertical angle is about 
10 degrees, the condition of a thin slice is realized. There- 
fore, the actual value should be a little above the inter- 
eats (Hp@)x, and is given to a good approximation by 

q. (3). 


AND CRANE 


values quoted by us, these corrections have been 
included. 
Scattering in the gas 

The directions of the incident and emergent 
tracks are obtained by fitting portions of circles 
to them; therefore, the direction obtained 
represents an average direction along its length, 
which is about 7 cm. The electron actually 
suffers small angle multiple scattering along its 
path in the gas, so that, if there were no scat- 
tering at all in the solid scatterer, we would still 
observe a small angular spread, due to the 
method of measurement and the scattering in 
the gas. The latter angles are added in random 
direction to those produced by the solid scatterer, 
so the measured result is the square root of the 
sum of the squares of the average angles pro- 
duced by the solid scatterer and by the gas. 
In fact, the gas can simply be considered as an 
addition to the thickness of the solid lamina. 
This correction has been calculated, and is small 
enough to be neglected, for all the scatterers 
used in the present work. For very thin scat- 
terers, such as the 0.0025-cm aluminum used by 
Slawsky and Crane,'® the correction is important, 
as Williams has pointed out. 


Error in measurement of energy 


The error in the measurement of the curvature 
of a track is due mainly to scattering in the gas, 
which makes it depart from circularity, or 
changes its apparent radius of curvature. 
Williams has derived a formula for this effect, 
and for tracks 7 cm long, in air at NTP, it 
reduces to p/p,=48 H, where p is the radius of 
curvature, due to the magnetic field alone, and 
p, is that due to scattering by the gas alone. 
The error for individual tracks in our measure- 
ments, where /7 ranges between 1260 and 2520 
gauss, is, according to the above, about 4 percent 
at 1260 gauss, and 2 percent at 2520 gauss. 
The average angle in the multiple scattering by 
the lamina varies inversely as the first power of 
the energy, so the error due to the above cause 
will remain symmetrical in the final result. 
Therefore, the estimated error, at 1260 gauss, 
is 4 \/n percent, where is the number of 


10 M. M. Slawsky and H. R. Crane, Phys. Rev. 56, 1203 
(1939). 
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tracks in the group (about 400). This is in the 
neighborhood of 0.2 percent. 


Energy loss inside the scatterer 


For each of our tracks, the arithmetic average 
of the ingoing and outgoing curvatures was 
used, since the energy difference was small, and 
one can, therefore, assume a linear decrease. 
This automatically takes care of the effect of 
energy loss, because, in the multiple scattering, 
the average angle depends linearly upon the 
energy. 


Path length in the scatterer 


Williams? and Rose" have called attention to 
the fact that, in scattering measurements, a 
correction should be made for path length in 
the scatterer. The path length” is approximately 


t(1+36), 


where ¢ is the actual foil thickness, and @ is the 
average angle of scattering. In terms of the 
average projected angle & this would become 


t(1+ 


if all angles concerned are small. This represents 
a minimum correction, since it assumes that all 
the tracks are incident upon the scatterer exactly 
perpendicularly. Actually, the tracks included 
in our data made angles up to 15 degrees with 
the normal as measured in the horizontal plane 
projection. The spread in angle in the vertical 
plane was much smaller, due to the geometry of 
the slit system through which the electrons 
passed. All facts considered, it seems fair to 
take ¢(1+24*) as the average path length. Thus 
the correction factor for the average projected 
angle of scattering is approximately (1+ 4). 


Personal errors in measurement 


If the errors in matching the circles to the 
tracks and measuring angles were random, or 
even symmetrical as to + and —, the net 
effect would be very small. Unfortunately there 
is no assurance that this is so, since the person 
measuring may have a tendency to make an 
error in a given direction. One can only guess at 


1M. E. Rose, Phys. Rev. 57, 556A (1940). 
2H. A. Bethe, M. E. Rose and L. P. Smith, Proc. Am. 
Phil. Soc. 78, 573 (1938). 


the amount of the error introduced. We believe 
that a generous estimate of the probable error 
in the final result due to such a systematic effect 


would be 3 percent. 


Calibration of instruments 


The errors involved in the calibration of the 
magnetic field, the size of the projected image, 
etc., combined, are expected to contain a prob- 
able error of about 2 percent. 


Combined effect of all errors 


In consideration of all the sources of error 
mentioned above we believe that a generous 
estimate of the probable error on our results is 
+5 percent. 


RESULTS 


Experimental results, with the values calcu- 
lated by means of the two theories, are given in 
Tables I and II. The theoretical values contain 
the necessary correction for the path length in 
the scattering lamina, the correction for scat- 
tering by the orbital electrons and the effect of 
cloud-chamber geometry (in the case of the 
Williams theory). Typical examples of plots of 
the data obtained together with Gaussian 
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Fic. 6. Scattering of electrons by Pb, 0.157 g/cm?. 


curves are shown in Figs. 3 to 7. Plots of all 
the curves obtained are not reproduced here 
because they would add little or nothing to the 
information contained in the accompanying 
tables. All the plots are Gaussian in shape, to 
within the experimental accuracy, as is to be 
expected. 

8 The experimentally observed @ was used in plotting 


the Gaussian. There would be a slight broadening of the 
curve if a theoretical value of @ were used instead. 
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384 OLESON, CHAO, AND CRANE 
TABLE I. Data on scattering by carbon (0.82 g/cm**). 
RATIO OF GovupsmMItT- RATIO OF 
NUMBER ENERGY OBSERVED SAUNDERSON OBSERVED 
OF INTERVAL OBSERVED OBSERVED WILLIAMS TO APPROXIMATE To 
TRACKS MEv (Hpa) yy, "1073 (Hpa) yy WILLIAMS (Hp9) G-S 
302 3-6 12.5° 2.27 2.56 89%, 2.51 01% 
243 6-8 8.4° 2.12 2.50 85% 2.45 87% 
297 9-11 7.0° 2.19 2.49 88°> 2.44 90% 
842 3-11 9.4° 2.20 2.51 2.46 
TABLE II. Scattering of 6- Mev electrons in various metals. 
RATIO OF GoupDsMIT- RATIO OF 
nd NUMBER OBSERVED SAUNDERSON OBSERVED 
__ OF OBSERVED OBSERVED WILLIAMS To APPROXIMATE To 
cM? TRACKS a (Hpa) Witttams (Hp) "1075 G-S 
Cu 0.174 543 9.0° 193 88% 93% 
Fe 0.173 435 9,2° 1.84 2.11 87% 1.98 93° 
Cd 0.111 633 9.4° 2.01 2.13 94% 2.01 100°% 
Pb 0.072 641 8.6° 1.72 2.07 83% 1.98 87% 
Pb 0.157 397 10.3° 2.68 3.21 84% 3.04 88% 
Pb 0.362 273 11.4° 4.03 5.12 79% 4.90 82% 
PbO 0.076 635 9.2° 1.76 2.08 85% 2.00 88% 
(0.070 Pb) 
Pb(C1O4)» 0.070 Pb 313 8.8° 2.21 2.75 81°¢ 2.62 84% 
Solution 
It should be noted that while most of the crystalline materials under nearly identical 


scattering materials listed in the tables are 
polycrystalline solids, two materials were used 
which were intended to be amorphous, PbO 
and Pb(CIO,)o..4 The purpose of this was to 
compare the scattering in crystalline and non- 


4 PbO in an amorphous form (massicot) was prepared 
in the following manner: lead hydroxide was first obtained 
by means of the reaction 


The lead hydroxide was then deposited on filter paper in 
several layers by forced filtering, care being taken to 
obtain uniform deposits. Several of these sheets of filter 
paper were carefully put together while wet and baked in 
the absence of air from 5 to 6 hours at 200 degrees centi- 
grade. This caused the formation of yellow lead oxide 
according to the reaction 


Pb(OH).+heat—~PbO+H.0. 


The resulting foil was used as the scatterer, appropriate 
corrections being made for the presence of oxygen and the 
filter paper. Lead perchlorate was prepared in the following 


manner: 
PbO+ 


The scatterer consisted of a cell having 0.2-mm glass walls 
containing a layer of lead perchlorate solution 2 mm thick 
with a concentration of 0.35 gram per cubic centimeter of 
lead, thus having a concentration of lead per square 
centimeter approximately the same as the thinnest lead 
foil used. The scattering power of all the atoms of the cell 
(including the glass walls) was taken into account, but 
was small compared to that of the lead. A preliminary 
report on the liquid scatterer was given at the Pittsburgh 
(1940) meeting of the American Physical Society: N. L. 
Crono)" kK. T. Chao and H. R. Crane, Phys. Rev. 58, 201A 
1940). 


conditions. The lead perchlorate solution was 
considered to be a particularly good approxi- 
mation to an amorphous scatterer, because there 
could be little order in the arrangement of the 
lead atoms in the solution and one could consider 
the scatterer to be very much like a “‘lead gas.” 


CONCLUSIONS 


The results seem to indicate that the scattering 
in carbon and other elements at the lower end 
of the atomic number scale is somewhat below 
the theoretically predicted values. These results 
are in agreement with our previous work on 
carbon,'® but diverge considerably from the 
values reported by Sheppard and Fowler,'® who 
found the scattering of 10-Mev electrons to be 
114 percent of theory (Williams) for carbon 
foils of 0.21 g/cm? and 0.61 g cm’. 

The experimental values for lead are also 
appreciably less than the theoretically predicted 
values, but the agreement with theory appears 
to be better than obtained by Sheppard and 
Fowler,'® who found 55 percent of theory 


%& N. L. Oleson, K. T. Chao, J. Halpern and H. R. Crane, 
Phys. Rev. 56, 482, 1171 (1939). 

°C. W. Sheppard and W. A. Fowler, Phys. Rev. 57, 
283 (1940). 
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(Williams) for a foil of 0.157 g/cm? and 71 
percent for a foil of 0.362 g/cm*. The lead 
perchlorate solution gives approximately the 
same results as the polycrystalline lead foils. 

The best agreement between theory and 
experiment is obtained for the case of the 
cadmium scatterer, an element in the middle of 
the periodic table. 

In conclusion, we believe that it is possible to 
state that there is evidence that the scattering 
for light and heavy elements is appreciably less 
than the theoretical values, while the middle 
elements appear to be in good agreement with 
theory. Although we have no theoretical pro- 
posals to make concerning these discrepancies, 
it might seem reasonable to expect the best 
agreement for the middle elements, because here 
the number of orbital electrons is large enough 
to justify the use of a statistical treatment of 
their distribution and also the atomic number 
is not so large that there is doubt as to the 
validity of the Born approximation. 

Our attempts to detect an effect of crystal 
structure gave negative results within our rather 
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Fic. 7. Scattering of electrons by a solution 
of lead perchlorate. 


broad limits of error. This, of course, does not 
preclude the possibility that a small effect exists. 
We wish to thank Professor Goudsmit for 
his continued interest and for the opportunity 
to discuss with him various aspects of the 
problem. This work was made possible by a 
grant from the Horace H. Rackham Fund. 
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Theory of the Magnetron. I 
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A complete calculation of space charge and field repartition is given for a magnetron working 
under steady conditions. Electrons leaving the filament gradually acquire an angular velocity, 
and for distances greater than a certain length L, these electrons describe spirals around the 
filament. This very important length L is defined by 

L? = —el/mwy'. 

7=current per unit of length of the filament, w7 =Larmor’s angular velocity. Under critical 
conditions, that is, when the magnetic field is just high enough to cut the anodic current J, the 
electron cloud rotates about the filament almost as a solid body with an angular velocity wy. A 
study of small oscillations with cylindrical symmetry shows that these oscillations have a proper 
frequency V2wz, and that the magnetron is able to yield an internal negative resistance for 
certain frequency bands near \2wy; this explains how a magnetron with one cylindrical anode 
can sustain continuous oscillations in an electric circuit. 


1. INTRODUCTION. GENERAL OBSERVATIONS ON 
THE PART PLAYED BY THE MAGNETIC 
FIELD—LARMOR’s THEOREM 


LLOWING Hull’s original work on magne- 


trons, a large number of theories have been 


formulated ; the majority, however, appear inade- 
quate or inexact. Certain authors considered the 
electronic motion while neglecting the space 
charge; others used the space charge computed 
by Langmuir for a diode without a magnetic 
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field. Such studies necessarily are highly inaccu- 
rate, since a powerful magnetic field considerably 
modifies the electronic paths and, consequently, 
the space charge. Contrariwise, authors who have 
endeavored to evaluate the space charge in t 
presence of a magnetic field do not appear to have 
provided practicable solutions. It therefore ap- 
pears necessary to consider the problem afresh 
and to compute directly both the space charge 
and the potential distributions in the magnetron ; 
without these essential data, the development of 
a coherent theory is impossible. 

As a preliminary study to the theory of the 
magnetron it is interesting to state some general 
results about electronic paths in an electrical 
field repartition with central symmetry, when a 
magnetic field is acting, on the system. Let us 
assume an electrical device where the potential 
V(r) is a function of the distance r to the z axis, 
while the magnetic field // lies parallel to the 
same axis. The equations of motion follow : 


al aVx 
me = —e— + poe, J] = —e— 
Ox or r 
(1) 
av y 
mij = —e— — woev, = —e— perl, 
orr 


where yo is the magnetic permeability in vacuum. 
By multiplying the first equation by —y and the 
second by x, adding and integrating, the integral 
of the moment of momentum is obtained: 


=m(xy— ye) + =c". (2) 


Let us introduce cylindrical coordinates 7, 6, 2 
around the z axis 


O/m=r6+ Hr =C; 
hence 
6=an+C/r’, (3) 


where 
wy = — (4) 


represents Larmor’s angular velocity. To de- 
termine the orders of magnitude, if H be meas- 
ured in gauss while all the equations are written 
in €.S.C.g.S. units, it is necessary to take: 


po=c, 


c=velocity of light; on the other hand, —e/m 


Y 


9 
x 


Fic. 1. 


= +5.3X10" e.s.c.g.s. since e is negative; then 
Wi =().884 x  — (5) 


which, for a field of 500 gauss, gives an angular 
velocity of 4.42 10°. 

If we revert to the equation of electron move- 
ment in polar coordinates r and 6, for 6, the 
simple result of (3) has been derived; for r, 


mi= — + mré", (6) 
r 


where the electric, Lorentz, and centrifugal forces 
will be recognized. If we replace woe/Z by the 
equivalent expression, — 2mw,, it is found that 


e ol’ 
j= (7) 
m or 


a result which may be expressed in this way: 


e oV 
-= —— —+(6-—wy)*?— wy’. (8) 
r my or 


Comparing Eqs. (3) and (8) we notice that the 
angular velocity comes in both formulae with the 
difference 6—wy which justifies the following 
treatment: In the xy plane, let us use a rotating 
axis OA, the angular velocity of which will be 
equal to wy; then the angle » between OP and 
OA as shown in Fig. 1(a) is 


n=0—wul, 
(9) 


and Eqs. (3) and (8) yield 


C e oV 
(10) 
ror mr Or 


The second equation contains only 7, which 
means that two opposite values, +C, giving 
opposite angular velocities +7 for each value of r 


A 


Ii 
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will yield similar orbits in the 7, 7 plane. The 
general solution is then represented by one orbit, 
rotated around the center O, Fig. 1(b), with 
Larmor’s precession wy, and by allowing the 
electron to run on this orbit (in the rotating 7, 7 
plane) with opposite velocity. This is the most 
general statement for Larmor’s theorem. The 
usual case considered in atomic structures refers 
to problems in which the angular velocity 7 on 
the orbit is very much greater than Larmor’s wy, 
so that the w,”? term in Eq. (10) can be neglected ; 
the rotating orbits are then similar to the 
unperturbed orbit obtained when the magnetic 
field is zero. 

In the problems connected with the theory of 
the magnetron, such an approximation will not 
be allowed, and it will prove necessary to take 
care of the w,* term. Let us, for instance, suppose 
that we want to study the motions of electrons 
inside of a cloud of constant density, p. The 
potential distribution V is given by 


lasoaVv 

AV=- -(—) =—4rpes.c.gs., (11) 
or 

V=—rZor’, 


where V and @1/0r are zero for r=0. Equations 
(10) vield 


n=C/r*, (12) 


The first condition gives the angular velocity as a 
function of r on each orbit, but different orbits in 
the cloud may correspond to different values of C, 
according to the initial velocity distributions for 
each special problem. 

For instance, we may suppose circular orbits, 
each of them corresponding to a certain C value. 
Such circular orbits will keep the electronic 
density p of the cloud constant; # being zero, the 
angular velocity is given by 


m 


e 


m 


(13) 


In each of the two + solutions the angular 
velocity is constant, independent of 7, and the 


electronic cloud is rotating as a solid body 
around the Oz axis. 

A constant linear velocity along the z axis may 
be added without any change in the preceding 
solutions, thus giving a possible type of electronic 
motions in the form of ‘“‘magneto-cathodic rays."’ 

At the limit p=0 one finds the well-known 
result that free electrons in a magnetic field may 
either stay at rest (6=0, sign —) or describe 
circular orbits with an angular velocity twice the 
Larmor rotation (6=2wy, sign +). 

Stable electronic clouds in magnetic fields are 
possible when their density p is smaller than a 
certain maximum value 


0 ps Pmy Pm= (14) 


If a cloud of higher uniform density had been 
artificially created, it would expand according to 
Eq. (12) (#>0), until it reached the limit p,.. 
Clouds of density smaller than p,, may be stable 
or expand or contract to the limit p,, the type of 
solution depending on the initial conditions of the 
cloud. 

It will be found that electronic clouds of such 
types play a very important role in the theory of 
the magnetron. 


2. CYLINDRICAL MAGNETRON—STATIC CASE 


The magnetron is assumed to consist of a 
filament of radius a and a cylindrical anode of 
radius b, the magnetic field 7 accurately paral- 
leling the filament (axis Oz). It is assumed that 
electrons without appreciable speed are emitted 
from the filament, and that the electric field on 
the filament is zero, provided the anode current is 
below saturation; hypotheses which have both 
been generally accepted since they were formu- 
lated by Langmuir. Distribution of the cylindrical 
space charge can thus be obtained ; the potential 
V is solely a function of r. The equations of 
motion are Eqs. (1), (2), (3), and (4) of the 
preceding section. 

The static case, which will first be considered, 
is characterized by the fact that V does not 
depend on time. 

The constant C of (3) is determined by the fact 
that electrons are emitted without speed from the 
filament, so that 6 is zero for r=a: 


6=wy(1—a?*/r’). (15) 


) 
) 
h 
r 
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Vo (b) 
<4 “1G. 2. Critica 


(4) 


The Lorentz force, due to the magnetic field /7, 
does not influence the velocity since it is perpen- 
dicular to the latter. Energy conservation ac- 
cordingly may be written: 


and this new constant is zero if 


V(a) =0 on the filament; (16) 


hence 


P+ (17) 


The equations of movement, therefore, may be - 


integrated directly from (15) and (17) without 
determining the potential distribution as a func- 
tion of r. 

In (17) one result is at once evident: since the 
velocity of rotation 6 is determined by the 
magnetic field, the kinetic energy at the distance 
» cannot be less than: 


2 2 


Hence, if the potential V(r) does not suffice to 
give the electron greater kinetic energy, the 
radial velocity * is annulled and the current is 
interrupted. Consequently 


r=0, —(2e/m) (18) 


This is Hull’s limiting value of the potential Vy at 
the distance 7, where the anode current is just cut 
off by the magnetic field. This limiting value 
Vo(r) is thus defined without the necessity of 
determining the distribution of the space charge 
or the potentials between the anode and cathode 
—a remarkable fact which discloses the possi- 
bility of finding the correct value of the critical 
potential Vo(r) without considering the space 
charge. 

Let us now come back to Eq. (7) in polar 


coordinates. Making use of Eq. (15) we find 


mi = — 6) 
= —edV/dr — 
= —(0/dr)[eV+ (19) 


The radial acceleration * is accordingly governed 
by an apparent potential function P(r): 


mi'= —e(0P/dr), 
(20) 


In the function P(r) a constant —mwy*a* has 
been added in order to reduce P(a) to zero on the 
filament r=a; thus the same function is made to 
appear as in the energy equation (17). A fact 
worthy of emphasis is the following: 

With the apparent potential P(r), the radial 
movement of the electron may be studied by 
means of Eq. (20) without considering the 
rotation w around the filament. These general 
results are valid in all cases either with or 
without space charge. 


3. STATIC SPACE CHARGE: CRITICAL POTENTIAL 


To obtain a potential distribution in a static 
state, it is necessary to introduce the space 
charge p(*), a function of the radius ry and inde- 
pendent of time. From the cylindrical symmetry 
of the svstem we have: 


AV = —4arp, IT=2nrpv,. (21) 
r or 


The current J, per unit length of filament, is a 
constant independent of r; v7, is taken from (17) 
and the following is obtained: 


as 2] 
r—)=-— 
or\ oar v; 
=— 


The conditions are the following: 


V(a)=0 on the filament. 
(0V/dr),..=0 no saturation. 


Since the non-linear Eq. (22) would require 
detailed discussion, it is preferable to start with 


( 
( 
1 
a 
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the simpler cases where the current I is sero. Two 
methods are available, according to Eq. (21): 


(A) 


which yields 


p#0, v,=0, 


m a*\? 
Vo(r) = P(r)=0. (23) 


e r 


Expression (23) completely satisfies the two 
limiting conditions of (22). 


(B) p=0, 


hence, 


B logr+C, (24) 


which corresponds to a definitely electrostatic 
potential without space charge. 

Let us now examine the distributions of the 
charges and potential in a magnetron when the 
potential of the anode (r=6) is rather low so 
that no current flows. 

The limiting case corresponding to the maxi- 
mum possible potential on the plate, with no 
current flowing, is obtained by taking 


V(b) = —(m/2e)wu?(b — a/b)? (18) 


as shown in preceding section, Eq. (18); it is the 
critical potential. The density is not zero, and 
space charges are present throughout the medium 
between the filament and anode; Eq. (21) gives 
their value: 


= (25) 


The radial velocity is zero throughout; the 
electrons follow circular trajectories, centered on 
the filament. The Lorentz and centrifugal forces 
are in exact equilibrium with the electrostatic 
force in Eq. (19), as is also evident from (20), 
since the distribution (18), (25) completely 
annuls P(r). 

In the case of a potential Vo(b) below the 
critical potential, there is obtained a charge 
distribution (25) extending from the filament 
(r=a) up to a certain cylinder r=d’. This par- 
ticular form of distribution then stops abruptly ; 
the potential distribution continues in the form of 
a logarithmic field at the cylinder 5’, and the two 
constants B and C may be determined by the 
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conditions: 


m a? 2 
——wy"{ } =B logh’+C 
( b’ 


(26) 


When the potential V(b) of the anode is pro- 
gressively lowered, the space charge is corre- 
spondingly restricted to a decreasing cylinder }’ 
around the filament; and, when V(b) is zero, the 
space charge disappears. 

Figure 2 illustrates the case of critical potential 
V(b); the different curves represent the space 
charge density p, the electrostatic potential V(r) 
and the apparent radial potential P(r). 

Figure 3 shows the potential distribution, V(r), 
for anode potentials lower than the critical value. 
The apparent potential energy of the electron is 
eP or —eP, ¢ designating the absolute value of an 
electron charge. Hence, 


P=V—-—Vo 


according to (20). The behavior of this function 
is easily described and is important to recognize 
inasmuch as the apparent energy P governs the 
electronic radial movements. 

In the critical state (Fig. 2), —P is identically 
zero. The electrons of the cloud forming the space 
charge have no radial velocity. 

When the magnetron is below its critical state 
(Fig. 3), the apparent potential —P rises between 
b’ and 6 in front of the anode. The critical 
distribution shown in Fig. 1 seems to be that 
described by Hull! in a brief note, for which he 
claims to have found good experimental evidence 
from space charge density measurements. 


Fic. 4. Values from 
Langmuir’s solution for 
the diode. 


Fic. 3. Potential dis- 
tribution, V(r) for anode 
potentials lower than the 
critical value. 


1A. W. Hull, Phys. Rev. 23, 112A (1924). 
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4. DIRECT-CURRENT CONDITION WITHOUT 
SATURATION 


Since the fundamental equation (22) has just 
been considered for the various cases where the 
current is nil, it is now desirable to determine the 
direct-current conditions. Equation (22) will be 
transcribed, with the unknown function taken 
as the apparent potential P(r) defined in (20) 
instead of the electrostatic potential V(r). 


or or 


If we replace Vo by its value (20) we have 


2mwn" a‘ 0 
or or 
2m\} 
= -1(-—") (27) 
e 


This equation is rigorous. The solution P(r) 
must be such that conditions (22) will be satisfied 
on the filament. 


P(a)=0, @P/dr=0, r=a 
on the filament (27). 
The solution of (27) can only be derived 


approximately; two distinct regions must be 
considered as limiting cases: 


A. Proximity to filament 

The second member remains constant and, on 
the filament, P is nil; the term in brackets must, 
therefore, be infinite. The required solution is of 


the type: 
P=A(r—a)", 


which gives 


2mwy? 
2a+nA(r—a)""! 
e 


2m 


+aAn(n— = -1( -— 


If we take n=4/3, the two first terms are zero; in 
view of the factor (r—a)"/? the last term remains 
finite and it is found that: 


P= Cer y—a<l. (28) 


This solution is valid only in the immediate 
vicinity of the cathode, while r—a remains very 
small. At increasing distances from the cathode, 
account must be taken of the fact that 0P/dr is 
no longer zero; and the magnetic term in wy 
must be dealt with. Let us first determine what 
happens at a short distance from the filament 
where the magnetic term still remains very 
small; the problem of the ordinary diode is 
re-encountered. 

The diode without a magnetic field has been 
dealt with by Langmuir ;? the potentials V and P 
are identical and Eq. (27) reduces to 


0 OP 2m\! 
-1(-—*) (29) 


For r>a, far from the filament, a solution is 
found in r} and the complete solution may be 
written: 

(30) 
where 8° is a function of r/a, calculated by 
Langmuir: 
r/ja_ 1 1.25 1.5 1.75 2 2.5 5 


0.045 0.116 0.2 0.275 0405 0.512 
4 5 6 7 8 9 10 15 oe 


0.665 0.775 0.818 0.867 0.902 0.925 0.94 0.978 1 


Equation (30) presents the same structure as 
(28) and correspondence is thus established in the 
vicinity of the filament. The potential P is 
obtained by replacing by (r—a)?/a. With 
close to a: 


= (r—a)*/ar. (31) 


Comparison with (30) shows that this ap- 
proximation is valid only up to r=1.25a; that is, 
in the immediate vicinity of the filament. This 
indication gives the limit of validity of (28), 
which should be replaced by Langmuir’s Eq. (30) 
when distances beyond the immediate vicinity of 
the cathode are involved. 

Figure 4 shows the shapes of the various curves. 


21. Langmuir, Phys. Rev. 2, 458 (1913). An approximate 
solution may be found: 


P=}(91)28 (- m 


it vields correct results for very small or very large r, the 
error is about 20 percent near r=4a 


( 
I 
t 
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B. Distant from filament 


At a considerable distance from the filament 
the Langmuir equation gives low values for the 
term (0/dr)[r(@P/dr) ]. Contrariwise, the mag- 
netic term within the brackets of (27) increases 
indefinitely. A moment therefore arrives when 
this first term becomes preponderant. At the 
limit, the second term may be neglected and the 


first term kept. An asymptotic solution P., valid . 


at a great distance, is: 


r>d«a, 


-2 
= +2) — /?—____ . (32) 


It is now necessary to discuss how to join these 
two extreme cases with a view to ascertaining the 
corresponding values of 7 and the form of 
function P(r) obtained. 

By comparing the equations, an expression 
—eI/mwy* appears, homogeneous to a certain 
length Z squared. The following gives an indi- 
cation of the order of magnitude: 


—e/m=5.3 X10" (e.s.c.g.s.), 
I=3X10°*J (J in milliamperes; J in e.s.c.g.s. 
units). 


If we calculate wy, according to (5) and express // 
in gauss: 


L?= —el/mwy® (33) 
The role played by the characteristic length L 
is very important and has now to be discussed. 
First, the condition stated at the beginning of 
Section 3 must be satisfied in order to justify use 


of the solution (32); this necessitates: 


By neglecting a, this condition reduces to 
r>L. 


For weak fields, which scarcely disturb the 
diode, the length Z is much greater than the 
dimensions of the bulb; if the field increases, 
however, the length Z decreases considerably to 
the order of magnitude of the interior dimension 
of the bulb and may even drop to very low values; 


2m 
<K——wu 


~)- “ mwyir® é 


Ve) 


Fic. 6. Magnetron 
characteristic. 


Fic. 5. Value of the 
apparent potential P,, for 
various values of L. 


for some hundreds of gauss, L is of the order of 
some tenths of a mm. 

The length Z then gives the magnitude of the 
distance where the Langmuir solution P,, (always 
valid in the immediate vicinity of the filament) 
has to rejoin with the solution P, which is valid 
at a great distance. 

As a rough criterion of the distance of coinci- 
dence, we take the condition: 


=P.(r), 
then 


3 (918?) (34) 


— [*e/2mwy'r? ; 


hence 
= L'. 


The curve representing the apparent potential 
P(r) is consequently known by the expressions 
P, (for r&L) and P,, (for r>L); it can, in any 
case, be plotted approximately. Figure 5 shows 
curves for various values of L between the 
cathode (r=a) and the anode (r=d). 

If the anode radius } is appreciably greater 
than the length L, the formula P,(R) of (32) 
gives the anode voltage, calculated from the 
critical voltage Vo(R); if, however, the length L 
is very great, the magnetic field being weak, 
Langmuir’s curve P, is used. Figure 6, ac- 
cordingly, shows the magnetron characteristic ; 
the internal resistance R is zero in the vicinity of 
the critical point. 

Reverting to the curves of Fig. 5, it should be 
noted that the electrons, after leaving the fila- 
ment (r =a), first reach the Langmuir region and 
then travel toward r=ZL through a region of low 
apparent potential energy, thereafter rebounding 
(approaching the anode) toward higher potential 
energies only slightly below the potential energy 
of the cathode. The expression — P(r) represents, 
but for the factor e= |e!, the apparent potential 
energy of the electrons, controlling their radial 


391 
U 
2 
i 


392 


movements (Eq. (20)). The speed of rotation 
around the filament is defined by (15). The 
curves —P(r) show that the electrons are ac- 
celerated from a to LZ and retarded from L to 0. 
The space charge thus is increased in the second 
region; close to the filament, Langmuir’s space 
charge applies, and, at a great distance, the con- 
stant space charge of (25) applies approximately 


r>L, (35) 


p= / 


This space charge density far away from the 
filament is just the density we have found in Eq. 
(14) as being the greatest possible density of 
electrons in a given magnetic field. A more 
detailed study would require calculation of the 
junction of the curves in the region L. The 
following, however, is pertinent: 

The electronic paths start radially from the 
filament; they are slightly curved in the region 
r=L and then turn around the filament for r>ZL, 
as shown in Fig. 7. This can be understood 
readily from the following considerations: Let V’, 
and Ve be the radial and rotational components 
of electronic velocity ; Eqs. (15) and (21) vield: 


Instead of computing the ratio J%/1’,, let us 
calculate 


po 


where po is the space charge density (25) in the 
region 7>L. If the radius a of the cathode is 
neglected, as a first approximation, it is seen that 
the ratio pV,/poVe equals —1 when r=L, which 
proves that V, and Vs are of the same order of 
magnitude at a distance L from the filament. It 
appears from this discussion that the junction 
between the two asymptotic solutions A (near 
the filament) and B (far away from it) has to 
take place about the distance L from the fila- 
ment. Series of approximations starting from 
both sides can be obtained and have already been 
published in a former paper.* 


pV, 


3L. Brillouin, ‘‘Theory of the magnetron,” Elec. 


Commun. 20, No. 1 (1941). 
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Fic. 8. Oscillator circuit. 


Electronic paths 
from filament to anode. 


Fic. 7. 


5. OSCILLATIONS IN A MAGNETRON WITH CYLIN- 
DRICAL ANODE: SOME GENERAL REMARKS 


Different types of magnetrons are practically 
used for generating electromagnetic oscillations ; 
some of them possess one anode, of cylindrical 
shape. Oscillations of large amplitude are often 
obtained from these devices, but large oscillations 
are rather difficult to discuss theoretically and 
the calculations will be given only for the case 
of small oscillations, which do not introduce too 
great a perturbation in the original charge 
density distribution. 

The theoretical discussion will be conducted 
with the following approximations: The formula 
giving the electric current has to be completed 
with an additional term, taking account of 
Maxwell’s displacement current. In more general 
problems, one should also take account of the 
magnetic fields induced by the current flowing 
inside the magnetron, but we shall assume the 
cylindrical symmetry to be obeyed for all 
quantities, in which case there is no induced 
magnetic field for a magnetron of infinite length. 
Magnetic fields would appear only outside of a 
magnetron of finite length and result in radiation 
of electromagnetic waves from the magnetron. 
This is, of course, a very small and secondary 
effect, causing some damping to the oscillations 
in the magnetron. 

Equations of Sections 1 and 2 will thus be 
modified in the following way: The charge 
density p and the radial electric field E will be 
considered as functions of r and ¢, an assumption 
which secures cylindrical symmetry at any time. 
Equation (11) becomes 


1a 
~—(rE)=4rp, E= (37) 


r or r 


|| 
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The current J, including displacement current, 
will be written 


1 dE 
pt — ~) =2mrpri+-—. (38) 
dt 2 dat 


The dielectric constant €9 in vacuum, has been 
taken 1 in e.s.c.g.s. units. As is well known from 
general proofs, introduction of the displacement 
term secures the constancy of J along the whole 
circuit, which means here that J is independent 
of the distance 7, and a function only of ¢. All 
other formulae from Section 2 still apply, and 
especially Eqs. (15), (19), and (20). Time 
derivatives may be taken either at a given place, 
or following the motion of an electron ; these two 
definitions will be distinguished this way : 


0/dt derivative at r con- 
stant or ordinary 
partial derivative, 
(39) 
d/dt=0/dt+v,0/dr following the motion 
of an electron mov- 
ing with the velocity 
Vr, Uo. 


All quantities being independent of 6, because 
of cylindrical symmetry, a term in 10/06 has 
been omitted in formula (39). Equation of 
motion (20), for instance, has to be written 


/dt? = —edP/dr. (40) 


A well-known integration method, which has 
proved very useful in most electronic problems, 
was omitted on purpose in the preceding sections 
for the reason that it leads to some difficulties 
which will be discussed now. The method con- 
sists in taking the derivative d/dt of the product 
rE, following the motion of an electron: 


d 0 
—(rE) =—(rE)+i—(rE) 
dt ot or 


dE 
(41) 
t 


and in using Eqs. (37) and (38). As noticed 
before, J(t) is independent of r, which allows easy 
integration : 


rE(r, t)=2 f I(t)dt. (42) 


to is the time at which one given electron has left 
the filament, and ¢ the instant when it reaches the 
distance r. 

Let us now show the difficulties involved in the 
use of this formula, by applying it to the static 
case. E is independent of ¢; hence 


E(r) =2I7/r, (43) 


where r=/—(y is the transit time for an electron 
going from the filament to the distance r. Sup- 
pose the magnetic field to be such as to cut the 
anodic current, according to condition (18); J is 
naught, but electrons in this case are moving in 
circular orbits about the filament, with no radial 
velocity at all; hence the transit time r is infinite, 
while E keeps a finite value. Many theoreticians 
attempting to use (43) without a former knowl- 
edge of the solution were led to entirely wrong 
conclusions. 


6. SMALL OSCILLATIONS IN A MAGNETRON WITH 
CYLINDRICAL ANODE 


After the magnetron has been studied in its 
static condition, without oscillations, all we need 
do now is to superimpose small oscillations on the 
static quantities. The amplitude of the oscilla- 
tions will be specified by a factor €; static con- 


_ tinuous quantities are characterized by a sub- 


script c, while alternating quantities will be given 
a subscript a. 


(V=V.(r)+eV,(r, t) electric potential, 
electric current in- 

dependent of r, (44) 

position of the elec- 

tron at the time ¢. 


r-(t) +er,(t) 


Because of the supposed smallness of ¢, all terms 
in &, &--- will be neglected. Let us now rewrite 
the most important formulae, starting with (20): 


or or 


m or 


d*r, d*re e fOVy e aV, 
( ) (45) 
dt? dit? m 


Potentials V») and V, have to be taken at the 
point r-+er, where the electron is located; hence 
the expansion = V(r.) /Or]. 
Constant terms and oscillating terms are easily 
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separated, and yield the following relations: 


d*re Volre) 0 
a m or or 
d*r, e e 
—— ( ) (47) 
om or? or? m Or 
Treating Eq. (42) in the same way, one gets 
eV. 
or? or 
t 
=2f telat; 
to 
hence two relations are derived 
OV. 
re— = —21(t—ty), (43) 
or 
OV. OV, 
Va (48) 
or or? 


Equations (46) and (43) between static quantities 
have already been obtained in the preceding 
sections and bring nothing new;; this is the result 
of our hypothesis of infinitely small oscillations. 
If oscillations of greater amplitude should be 
considered, terms in é, &--+- could not be neg- 
lected any more, and would appear in the static 
formulae (corresponding to detection phenomena) 
thus perturbing the static regimen. 

Equations (47) and (48) rule the magnetronic 
oscillations and are linear with respect to the 
unknown quantities r, and V,. Eliminating V, 
one easily gets 


OV. 

-— —|-- I,dt, (49) 
dt® or? or, mr, 
where Vo and V, result from the study of the 
static case. 

In order to know how the magnetron is able to 
react on the oscillating circuit connected to it, we 
need to calculate the internal impedance of the 
magnetron; we thus want to find the relation 
between the alternating current J,(t) flowing 
through the magnetron, and the alternating 
potential V, on its anode. Let us suppose J,(¢) to 
be given, then Eq. (49) is a linear differential 
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equation with a right-hand term. As is the rule 
for such equations, we have to study first the 
linear differential equation without the second 
member; it is of the ordinary oscillating type, and 
indicates a proper frequency w(r) as a function 


of r. 
or? 


lav. 


(50) 
r or 


Let us study first the order of magnitudes of these 
internal proper frequencies. In a magnetron on 
critical conditions, the potential repartition is 
given by Eq. (23) 


Vo=Ve= (51) 


hence 


w*(r) = w(r)=Vlwy. (52) 


Such critical conditions, however, correspond to 
no anodic current (J-=0) and infinite transit 
time 7, and cannot be used for sustaining oscil- 
lations. We may, however, suppose the magnetron 
to be very near this critical state, the continuous 
current J, keeping a small value, so that the 
average potential V.(r) will differ only slightly 
from V(r). The proper frequencies w(r) will then 
keep values very near V2wy. 

Returning to Eq. (49) we notice its similarity 
to the equation of a harmonic oscillator acted 
upon by an external force represented by the 
right-hand term. A remarkable feature is that 
this equation contains no damping term, but this 
is only the result of the initial assumptions 
outlined in Section 5. A complete theory, in- 
cluding magnetic oscillating fields, should yield 
equations with a damping term due to outside 
radiation, as is immediately realized for physical 
reasons. 

We are thus allowed to proceed as though such 
a damping were actually found in the equations, 
and search only for the forced oscillations, 
neglecting free oscillations [solutions of Eq. (49) 
without right-hand term]. Let us thus suppose 
conditions very near the critical state and call 


=t—to (53) 


the transit time for an electron running from the 
filament (r=a) to the distance r. Supposing 
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Fic. 9. Approximate values of Ro. 


sinusoidal oscillations, we write down 


J,>0 real, 


Ja | Jo wr 
I,dt = 2— sin—e'+'¥, 
tw w 2 
¢=—wr7/2. 


Equation (49) takes the following form and 
defines, for each distance r., the amplitude r, of 
electronic oscillations. 


4eJ, 


—w*) = (55) 


Phase ¢ of these oscillations is directly related to 
the transit time r. 


From the r, values, one easily calculates V, 
using Eq. (47). 


OV. m 
or or? or? oe 


(56) 


V> and V, are very nearly equal to each other 
(almost critical state) and formula (56) ap- 
proximately yields 


Ore rLw?(r) —w? 
Va= {--- (58) 
r)—w 


One must, of course, suppose V, to be zero on the 
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filament r=a. These simpler formulae result 
from the assumption of a magnetron working 
very near its critical point. On a more general 
case, formulae would be of a more complicated 
type, but one very important point would remain 
unchanged, namely, the relation between phase 
angle ¢ and transit time r. 

Having obtained expressions for alternating 
current J, and potential V,, one immediately 
finds the internal impedance of the magnetron 


_(—--) *, (59) 
Je w*(r) —w* 


¢=—wr/2. 


The question is now to discuss this expression and 
to see under which conditions it may yield a 
negative real part, indicating a negative re- 
sistance term in the magnetron. It is well known 
that, in order to be able to sustain oscillations 
in an electric circuit with positive resistance, an 
electronic tube must play the role of a negative 
resistance. Depending on secondary effects in the 
electronic tube, one may connect this negative 
resistance in series with the electric circuit or in 
parallel with it. Magnetrons are actually used 
with the parallel connection, according to 
Fig. 8. 


7. THE MAGNETRON AS A SOURCE OF 
ELECTRICAL OSCILLATIONS 


Referring to formula (59) we notice, first, that 
the denominator is zero for w* equal to w*(r), 
which seems to indicate the possibility of infinite 
internal impedance R,. This, of course, is the 
result of some simplifying assumptions discussed 
in Section 5. A complete theory would include 
direct radiation from the moving electrons, thus 
giving a damping term s in the equation, and the 
impedance R, would have a denominator avoid- 
ing infinity like 


{[w*(r) —w? 
At any rate, large values of R, will only be 


obtained for 


w= w(r) Vlwy. (60) 


| | 
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Now we determine the real part of R, 


2w 
Re(R,) = sinwr. (61) 


w*(r) —w*\a? 


Figure 9 gives an approximative visualization of 
such a function. This will be negative in the 
following cases: 


w<w(r), sinwr <0, 
r<wr<2r, (62) 
3a <wr <A4n, 

and 
w >w(r), sinwr>O, 
0<wr<z, 
2a <wr<3rn, 


(63) 


We recall that transit time 7 in a magnetron may 
reach values much greater than the period of 
oscillations; the best values of 7 in both cases will 
correspond to sinw7 equal 1, that is, 


wr 


(64) 
(65) 


w’<w*(r), 
w*>w*(r), 


One thing may appear rather surprising in the 
preceding results, and that is the role played by 
frequencies of the order of V2wy. If we refer to the 
electronic motions in the static magnetron, we 
notice that the electrons move around the 
filament with angular velocities ranging from 0 
(near the filament) up to wy at a great distance 
from the filament [Eq. (15) ]. None of these 
frequencies happens to play a role in the opera- 
tion of sustaining oscillations, and the v2wy 
frequency is noticeably higher. This can be 
understood, as the angular velocities of the 
electrons cannot be observed in a magnetron 
whose anode is built as a complete cylinder 
around the filament. The only motions which 
may be noticed as giving a change in the anodic 
current are radial motions. The oscillations 
result, in this case, in electronic motions, where 
all electrons located on a certain cylindrical 
laver (at distance 7.) move together either to or 
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fro, the whole cylindrical layer expanding or 
narrowing at one time. This is the type of 
motion corresponding to the frequency V2wy, and 
it has nothing to do with the rotational fre- 
quency w. 

Another point to be emphasized is the great 
variety of conditions under which the magnetron 
is able to sustain oscillations. Conditions (62) and 
(64), for w<w(r), and (63) and (65) for w>w(r), 
offer a great number of different solutions, a 
characteristic almost never found in the dis- 
cussion of electron tubes. 

Where does the energy necessary for sustaining 
oscillations in the outer circuit come from? All 
our calculations refer to a magnetron very near 
its critical conditions (J=0). Such a magnetron 
stores up energy in kinetic form from the 
electrons rotating about the filament. In the 
process of sustaining oscillations in an outer 
circuit, some energy will be taken out of this 
reserve, thus slowing down the electrons; in 
order to keep the electron cloud in rotation, a 
small current J must be allowed to flow between 
anode and cathode, and the power JV supplied 
by the battery will restore the necessary energy 
and keep the cloud in steady rotation. 

Our conclusion, that such a magnetron should 
be operated near vV2wy, fits with the results 
obtained by Blewett and Ramo?‘ in a recent 
paper and is checked by experimental facts. 
Blewett and Ramo, however, failed to get the 
conditions for negative resistance, a difficulty 
they themselves emphasized at the end of their 
paper. 

The first part of this paper, including Sections 
1-4, was written in January, 1939, but its 
publication was delayed because of present cir- 
cumstances. Several copies, however, were circu- 
lated; the author has to thank M. M. Blewett and 
Ramo for their courtesy in making reference to 
this unpublished paper. 


‘J. P. Blewett and S. Ramo, Phys. Rev. 57, 635 (1940). 
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A detailed study has been made of the lattice imper- 
fections which are assumed to account for plastic de- 
formation. Calculations of the strain energy associated with 
these dislocations have been made in the following cases: a 
dislocation in a uniform shear stress, two dislocations in an 
infinite medium, and a dislocation near a surface. The force 
acting on a dislocation is found by taking the gradient of 
the strain energy. A force is found which tends to attract 
dislocations toward the surface of a specimen. It is shown 
that about twice as much energy is required to produce a 
certain amount of slip inside a solid as is required to produce 


the same amount of slip at the surface. The energy required 
to produce a dislocation is found to be several electron volts 
per atomic plane, the exact amount depending on where it is 
located and how it was produced. Finally the energy stored 
in a material during work hardening is calculated by 
assuming that the dislocations are arranged in a regular 
two-dimensional lattice in the material. The density of 
dislocations found for severely work hardened material 
agrees with the predictions of other investigators. Numerical 
values found for the energy stored during work hardening 
are in agreement with experiment. 


INTRODUCTION 


F a single crystal is subjected to a small 
tensile stress it will elongate according to 
Hooke’s law. If the component of the shear 
stress acting along a definite crystallographic 
direction exceeds a certain critical value the 
crystal will also undergo considerable plastic 
deformation. For most crystals the value of the 
critical shear stress is about 10’ dynes per cm?. 
The permanent deformation apparently results 
from a sliding of sections of the crystal over one 
another. The sections are bounded by certain 
close-packed crystallographic planes and the 
slip takes place along the direction mentioned 
above. This slip direction is usually a direction 
in which the atoms are closely packed. A micro- 
scopic examination of a single crystal which has 
been elongated by several percent shows clearly 
the step-like discontinuities produced at the 
surface when slip occurs. The sections bounded 
by slip planes are several microns thick at room 
temperatures. 

In a series of articles appearing in the Journal 
of Applied Physics, Seitz and Read! consider the 
various theories which have been used to discuss 
plastic deformation. They conclude that the 
theory of ‘“‘dislocations’’ as developed by Taylor* 
yields the most satisfactory description of plastic 
phenomena which has been given to date. The 


*Rackham Postdoctoral Fellow of the University of 
Michigan. 

1F. Seitz and T. A. Read, J. App. Phys. 12, 100, 170, 
470 (1941). 

?G. I. Taylor, Proc. Roy. Soc. 145, 362 (1934). 


present paper represents a further attempt to 
develop and to test the validity of dislocation 
theory. The predictions of dislocation theory 
which are considered in this paper are in good 
agreement with experiment. 

According to Taylor, slip in a solid results 
from the motion of a certain type of lattice im- 
perfection through the crystal. Taylor has used 
these ‘‘dislocations” to discuss qualitatively the 
hardening which occurs when crystals are 
plastically deformed. Recently Mott and Na- 
barro® using dislocation theory have discussed 
precipitation hardening. 

In order to make clear the nature of a disloca- 
tion we consider a line dislocation in a simple 
cubic crystal. In Fig. 1 we have indicated the 
positions of the atoms in a crystallographic 
plane which is normal to the slip plane. The 
axis of the dislocation is perpendicular to the 
plane of the figure. Thus parallel atomic planes 
above and below the plane of the figure will have 
the same distribution of atoms about the center, 
or axis, of the dislocation. If a dislocation moves 
along its slip plane across the entire specimen, 
then Taylor has shown that the material above 
the slip plane will be displaced relative to the 
material below the slip plane by an amount 
equal to the lattice constant of the material. 

An examination of Fig. 1 indicates that the 
distortion is just that which would be produced 
by inserting extra half-planes of atoms into 


3N. F. Mott and F. Nabarro, Proc. Phys. Soc. 52, 86 
(1940); F. Nabarro, Proc. Roy. Soc. 175, 519 (1940). 
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an originally perfect crystal above the slip 
plane. Following Taylor we shall call this a 
positive dislocation. If the additional half-planes 
are inserted below instead of above the slip plane 
we shall call the dislocation which results a 
negative dislocation. Thus, if we are dealing 
with a positive dislocation we shall find that the 
material above the slip plane is compressed while 
the material below the slip plane has been 
dilated. 

In the first section of this paper we shall 
calculate the force which acts on a dislocation 
when the material containing the dislocation is 
subjected to an external shearing force. Next, 
we shall consider two dislocations in an infinite 
medium. We shall determine the energy of the 
system and we shall evaluate the force acting 
between the dislocations. In the third section 
we shall discuss the case of a dislocation near a 
surface. Here too we shall calculate the energy 
of the system and the force acting on the dis- 
location. Finally, we shall evaluate the energy 
stored in a material when it is work hardened. 
All calculations made in this paper will use the 
theory of an elastic, isotropic continuum. 


DISLOCATION IN UNIFORM SHEAR STRESS 


\ According to dislocation theory an annealed 
specimen contains only a few line dislocations. 
If we apply a sufficiently large external shearing 
force to the specimen it will deform plastically. 
This deformation may take place either by the 
motion of existing dislocations or by the forma- 
tion and motion of dislocations. In either case 
the force produced on an isolated dislocation 
by the external shearing force would be of in- 
terest. We shall, therefore, consider a line dis- 
location in a material which is also subjected to 
a uniform two-dimensional shear stress. This 
uniform shear stress is produced by the external 
force. 

We shall first calculate the strain energy of the 
system. In order to calculate the strain energy 
we must know the stresses acting. It can be 
shown that the solution of certain two-dimen- 
sional problems of elasticity can be obtained by 
finding a suitable Airy stress function,‘ x. This 
stress function must satisfy the differential 


*S. Timoshenko, Theory of Elasticity (McGraw-Hill 
Book Co., Inc., 1934), p. 25. 
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equation 


dy"? 


dy* 


and must yield stresses or displacements which 
satisfy the boundary conditions of the problem. 
The stresses are given by the following relations: 
o,=0°x/dy"; Try= —0*x/dxdy. (1) 

The Airy stress function appropriate for this 
problem is® 


Fic. 1. Dislocation lying in a simple cubic crystal. 


x= —Bxy—Dy log(x*+*)*. (2) 


The first term represents the uniform shear 
stress; the second, as one can see from the dis- 
placements given below, represents a positive 
dislocation. The stresses resulting are 


o2= 

The stresses due to the dislocation die off 
inversely as the first power of the distance from 
the center of the dislocation. The displacements 
can now be calculated from Hooke’s law. They 
are 


(3) 


m—1 D «x 
mG 
(4) 
(x?-+y2)! Dx? 
y= x . 
2mG 2G x*+y" 


5 This dislocation differs from the one used by Taylor. 
Dislocations of this type were probably first used to discuss 

lastic deformation by J. M. Burgers and W. G. Burgers. 
}. Mt Burgers, Proc. K. Akad. Amst. 42, 263, 378 (1939). 

e shall discuss the differences between these two types 
of dislocations near the end of the next section. 


Ss = 


THEORY OF PLASTIC DEFORMATION 399 


Here u is the displacement in the x direction 
and v is the displacement in the y direction. G 
is the shear modulus of the material, and m is 
the reciprocal of Poisson’s ratio. If we measure 
the angle appearing in the equation for u from 
the plus x axis then the displacements are dis- 
continuous along this axis. The material just 
below the positive x axis is displaced in the x 
direction with respect to the material just above 
the axis by an amount 


\=2r(m—1)D/mG, D=mGd/2xr(m—1). (5) 


According to Eqs. (3) and (4) the stresses and 
the displacements are singular at the center of 
the dislocation. Physically this means that an 
elastic continuum theory is not capable of 
describing the situation near the center of a 
dislocation. We shall, therefore, in all problems 
surround the axis of each dislocation by a 
cvlinder of radius 7». We shall choose 7» so that 
the maximum strain at the surface of the cylinder 
is one-tenth. We shall consider only situations 
in which all dislocations, impurity atoms, 
precipitates and other singularities are separated 
from the dislocation in question by distances 
which are large compared with ro. For such 
systems the stress at ro will be, to a good ap- 
proximation, just the stress produced by the 
dislocation. Furthermore, the interaction ener- 
gies and the forces exerted by other singularities 
on the dislocation will depend upon the strain 
energy in the material outside 7o since the strain 
energy inside ro will not change appreciably if 
the dislocation is moved relative to the other 
singularities. 

Let us calculate ro in the case of copper. We 
shall take A of Eq. (5) to be the translation 
which would carry an atom of copper from a 
cube corner to the center of a cube face. This 
unit slip distance is 2.55X10-* cm in copper. 
Using G=4.53 X10" dynes per cm? and m = 2.941 
we find from Eq. (5) that D=2.78X10* dynes 
per cm. The shear strain is found to be the 
largest strain. Using Hooke’s law and substitut- 
ing from Eqs. (3) we find 


Yey = Tey/G= +Dx(x?— y*)/G(x*+y"*)?, 
where yz, is the shear strain. Putting in the 


limiting value of one-tenth for y,, and inserting 
the maximum value that the expression on the 


right has on the surface of the cylinder we find 
to=D/Y2y max G=6.15+10-* cm for copper. (6) 


The strain energy of a slab of material of 
unit thickness in the z direction is® 


w= f 


(7) 


where m is the reciprocal of Poisson's ratio and 
where we have assumed that the material is in 
a state of a plane strain. Equation (7) gives the 
strain energy in any orthogonal two-dimensional 
coordinate system. In this equation ds is the 
element of area and a, o2, and 712 are the stresses 
in the coordinate system chosen. Substituting 
from Eqs. (3) we find 


9 


1 D? 
f f axiy| 
2G x?+y? 


2BDx(x*?—y*) 
(xty2)2 


It is apparent from Eq. (8) that if either the 
uniform shear stress or the dislocation is present 
in the specimen it will contain strain energy. 
If both the uniform shear stress and the disloca- 
tion are present simultaneously, the strain energy 
is not simply the sum of the strain energy due 
to a uniform shear stress and the strain energy 
due to a dislocation. In addition to these self- 
energies we have a cross term, an interaction 
energy. The first two terms of the integral (8) 
are self-energy terms while the last term is the 
interaction energy in which we are interested. 
In order to get a sensible result one must be 
careful about the limits of integration. The inter- 
action energy can be written as 


2BD — x(x? — 
at > —_——dxd 
J 


et+€ 
+f JS 


(8) 


In this expression C is a large distance which is 
*E. Trefftz, Handbuch der Physik, Vol. VI, p. 79 (1928). 
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allowed to approach infinity after the integra- The force acts in the plus x direction on this 
tion has been carried out. The result then positive dislocation. 
becomes Using the value of D for a unit dislocation in 
Wint= —2BDE/G. (10) copper and assuming that the uniform shear 
stress is 10’ dynes per cm (i.e., the value of the 
critical shear stress for copper) we find the 
force to be 2.22 X10~-* dvne for a segment whose 
length is the distance between successive atomic 
planes which are perpendicular to the axis of the 
dislocation. This force is very small compared 
(11) with chemical forces. 


The distance & essentially measures the position 
of the dislocation with respect to the field of 
uniform shear. The force acting on the disloca- 
tion is given by the negative of the derivative of 
this expression with respect to &. 


F=+2BD/G. 


Two PARALLEL LINE DISLOCATIONS 


In this section we shall calculate the strain energy associated with two parallel line dislocations 
in an isotropic, continuous medium. Quantitative values of the strain energy and of the force acting 
between the dislocations will be given. Finally, Taylor’s treatment of dislocation theory will be 


discussed. 


The Airy stress function for two unlike dislocations of equal strength in an infinite medium is 
x =D ty log(x*+y*)!— (y— yo) (x yo)? }}}. (12) 


We have assumed that the negative dislocation has its center at the origin of our coordinate system. 
The positive dislocation has its center at xoyo. Each dislocation has associated with it a slip plane 
which passes through the center of the dislocation and is parallel to the xz plane (zs is normal to the 
plane of Fig. 1). The stresses obtained from (12) are 


(x?+-?)? ({x—x0}*+ {y—yo}*)? 

(x?+-?)? ({x—xo}?+ {y—yo}?*)? 


From Eqs. (12) and (13) one sees that if the two dislocations coincide they will neutralize one another 


leaving an unstrained crystal. 
Substituting the stresses given in Eqs. (13) into expression (7) we obtain the strain energy for an 
infinite slab of material which is of unit thickness in the z direction: 


=—| | dxd + 

2x (x? (x? — 9") (3 — x0}? + {y—yo}*) 


90) Bx? +97) (1x — — {y—yo}*) 


» (14) 


where 


r=(x*+y")!, 


I 
f 
I 


1S 
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Let us surround the axis of each dislocation by a cylinder of radius ro and proceed with the integra- 
tion of Eq. (14) over the portions of the slab which lie outside these cylinders. The terms in Eq. (14) 
can be divided into two classes. First, there are the self-energy terms which contain r; or 72. Second, 
there are the interaction terms which contain both 7; and re. The self-energy terms can be integrated 


immediately and yield 
Weert = log—, (15) 
mG ro 


where L is essentially the radius of our cylindrical specimen which has the two dislocations near its 
center. We introduce new Cartesian coordinates into the interaction terms taking the origin midway 
between the two dislocations. The new coordinates ¢ and 7 are given by 


x= cos¢—» sind+3R cosd, y= sing+n cos¢+}R sing, (16) 
R=(xo?+y0?)!, tand=yo/Xo. 


where 


The interaction terms can then be written: 
D* 4(m— ps 7? cos*@ — } R? 
cos26+ ¢4(n? cos2o— cos2) + — cos2¢ — n?R*4+2 cos2¢ ]+[3R‘cos2¢/16 }) 
cos26+ —44+2 cos2¢ |+n? R*L4+cos2¢ cos2¢/64}) 
In Eq. (17) we have omitted terms that are odd in ¢ or 7 since these would obviously not contribute 
to the integral. If we now integrate Eq. (17) over the region outside the cylinders we obtain the 


interaction energy. Adding this interaction energy to the self-energy we find that the total, elastic, 
strain energy of the positive and negative dislocation in an infinite medium is 


R 


= — 32 18 


(17) 


Expression (18) is approximate since we have slab whose thickness is the distance between 
neglected small terms of order (m*Gd?/4x(m—1)*) successive atomic planes which are perpendicular 
X(ro?/R*). The force acting on the positive to the axis of the dislocation. In Figs. 3 and 4 


dislocation at xovo has the components 
1 


—1) R’ 
(19) 


sin 2@ 
R 


Numerical values of the strain energy and the 
force acting between dislocations have been 
calculated for various distances of separation in 
copper. The results are given in Figs. 2, 3, and 4. 
In Fig. 2 we have given the strain energy in a 


we have given the force acting on a dislocation 
segment whose length is equal to the thickness 
of the slab just mentioned. In Fig. 4 we have 
given that component of the force between two 
dislocations which acts along the slip direction. 
We have assumed that the slip planes of the two 
dislocations being considered are separated by 
by 10-* cm. 

We are now in a position to discuss some of the 
differences between the dislocation theory de- 
veloped by Taylor and the dislocation theory 
used in the present paper. Imperfections of the 
type shown in Fig. 1 can be obtained in elastic 
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ELASTIC ENERGY (EV PER PLANE) 


DISTANCE OF SEPARATION (CM) 


Fic. 2. Strain energy for various distances 
of separation in copper. 


theory if one uses the following stress function: 
x =Ay log(x?+y?)!+ Bx tan“'(y/x). (20) 


The strength and nature of the dislocation repre- 
sented by the above function are determined by 
the values assigned to the constants A and B. 
If we put A=B=—)G_r then we are dealing 
with a positive Taylor dislocation whose slip 
distance is \. If we take A = —mGX,22r(m—1) 
and B=0 we are dealing with a positive Burgers 
dislocation whose slip distance is also X. 

We have not used Taylor dislocations for the 
following reason. Suppose we have a block which 
is finite in the x and y directions and which 
contains a dislocation. Let us take the surface 
of the block to be free from external forces. 
Then it can be shown that the boundary condi- 
tions which apply to this problem can only be 
satisfied if one takes the dislocation to be a 
Burgers dislocation. 

In any dislocation theory it is important to 
decide how a dislocation will move along its 
slip plane under the influence of various stress 
fields which may be present in the material. 
Taylor treats the problem in the following way. 
Consider a positive dislocation located at the 
origin of a coordinate system. The xz plane will 
be taken to be the slip plane and the x axis the 
slip direction. Taylor, Burgers, and Mott and 
Nabarro assume that this positive dislocation 
will move in the plus x direction if the total 
external shear stress acting at the origin in that 
direction is greater than zero. If the above as- 
sumption is valid it would be most reasonable 


10°. 


FORCE ON DISLOCATION 
(OYNES PER ATOMIC PLANE) 


DISTANCE OF SEPARATION (CM) 


Fic, 3. Force acting on dislocation segment. 


to assume further that the force tending to 
move a dislocation in a given direction is pro- 
portional to the external shear stress which acts 
in that direction at the center of the dislocation. 
This last assumption can be shown to be in 
error in the following way. If the assumption is 
used to calculate the force acting between two 
Burgers dislocations the result found does not 
agree in form with the force obtained by calcu- 
lating the strain energy of the system. If one 
considers Taylor dislocations in various stress 
fields then the assumption gives the correct 
form for the forces, but the constant of propor- 
tionality seems to depend upon the particular 
stress system being considered. For example, 
one finds that 

2F4,ta=Fp/ tp. (21) 


Here F, is the force acting on a positive Taylor 
dislocation in a uniform shear stress; 74 is the 
value of the uniform shear stress at the center of 
the dislocation. Fx is the force acting between a 
positive Taylor dislocation and an equal nega- 
tive Taylor dislocation; rz is the shear stress 
produced at the center of the positive disloca- 
tion by the nearby negative dislocation. These 
forces are obtained by calculating the strain 
energy of the system considered. 

In spite of the discussion which we have just 
given it can be seen that the results of Taylor’s 
theory are qualitatively correct. This is true 
because Taylor assumed an interaction between 
two dislocations which is essentially of the correct 
form (i.e. the force dies off inversely as the first 
power of the distance between the dislocations). 
The error resulting from the use of the wrong 
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+6x10° 


Fic. 4. Force be- 


tween two dislocations 
in direction of slip. 


(OYNES PER ATOMIC PLANE) 


FORCE ACTING IN SLIP DIRECTION 


| 


DISTANCE OF SEPARATION OF DISLOCATIONS (CM) 


type of dislocation causes Taylor to overestimate 
the critical shear stress by a factor of about 
4(m—1)/m. For copper this factor has the value 
2.64. Our discussion indicates that Taylor’s 
method of determining when a dislocation will 
move may also be in error. If this is true then 
Eq. (21) allows us to estimate the error made. 
Such an error would cause Taylor to under- 
estimate the critical shear stress by a factor of 
about two. Thus Taylor’s results are not in 
error by more than a factor three. 


A DISLOCATION IN A CYLINDER 


There are two ways in which dislocations 
could be formed. First, dislocation pairs’ could 
be created inside the specimen. Each dislocation 
pair would consist of a positive and a negative 
dislocation. Secondly, single dislocations may be 
formed at the surface of the material. In order 
to consider these possibilities in detail we shall 
calculate the strain energy and the force acting 
on a dislocation near a surface. Finally we shall 
discuss briefly the annealing process in metals. 

Let us consider a dislocation in a circular 
cylinder of radius p. We shall suppose that the 
dislocation is not in general at the center of the 
cylinder. We shall take the axis of the dislocation 
parallel to the axis of the cylinder. The sym- 
metry of this problem is that of two eccentric 
cylinders. Problems of this type can be con- 
veniently treated using bipolar coordinates. 
G. B. Jeffery’ has used this orthogonal system 

7G. B. Jeffery, Trans. Roy. Soc. 221, 265 (1920). See 


also E. Coker and L. Filon, Photoelasticity (Cambridge 
University Press, 1931). 


of coordinates to treat two-dimensional elastic 
problems and we shall use many of his results. 
The curvilinear coordinates are defined by the 
relation 


x+i(y+c) 
(22) 
x+i(y—c) 

where ¢ is a real positive constant and x and y 
are Cartesian coordinates. Solving for x and y 


we find 


atis=log| 


x=c sin8/(cosha—cos@), 
y=c sinha, (cosha—cos§). (23) 


The points 0,+c and 0,—c are the two poles of 
the system. It is clear that if we keep 8 fixed 
and change a by an amount da we shall traverse 


a distance 
da, h=cda, (cosha—cos8) 


similarly if we keep a@ fixed and change 8 by dg 
we shall move through the distance d@/h. 
Figure 5 shows the two orthogonal families of 
circles that form this curvilinear coordinate net. 

Let us suppose that the dislocation is at the 
pole which lies at 0,+c. The surface of the ma- 
terial will be represented by the circle on which 
sinha=c/p. If the surface of the cylinder is to 
be free from external forces the boundary condi- 
tions which the stresses must fulfill are 


o2=0, Tap =0 (24) 
on the cylinder having sinha=c/p. We shall try 
a stress function which we believe to be ap- 
propriate for this problem. We put 


D sing 
{a+} exp[—2(a—ay)]}, (25) 
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where a is defined by sinha;=c/p. If the first 
term of this expression is transformed into 
Cartesian coordinates we find that it represents a 
positive dislocation at 0,+c¢c and a negative, 
“image” dislocation at 0,—c. This stress function 


{2(cosha—cos8) exp[ —2(a— a) ]—sinha[1—exp[ 


KOEHLER 


is a linear combination of terms which Jeffery 
has shown satisfy Airy’s equation. Jeffery has 
also given relations which enable one to calcu- 
late the stresses in bipolar coordinates once the 
stress function is known. Using these we find 


(26) 


D sinha sing 
D sing 
+ 
c 
D cosg 
tas=— (cosha — cos@) {1 —exp[—2(a—a1) ]}. 


c 


It is at once evident that these stresses satisfy the boundary conditions given above. 
We can now calculate the strain energy. Substituting the stresses into Eq. (7) we obtain 


D? 


2G 


f f 2(m —1) exp[ —4(a— a1) ]+m cos*8 {1 —exp[ —2(a—a,) J}? 


2(m—2) sinha exp[ —2(a—a1) ]{1—exp[—2(a— a) ]} 


cosha — cos8 


4 (m—2) sinh?a {1—exp[—2(a—a,) ]} 
(cosha —cos8)* 


In this equation we do not integrate @ to plus infinity because, just as in our previous problems, we 
find that elastic theory is not able to describe the distortion of the material near the center of the 
dislocation. Instead we put sinhag=c/ro. Carrying out the integrations we find 


mG)? 
W =———_ 
4x(m—1)? 


—1)a+3(m-—2) sinha {sinha —cosha} +( 


3m—4 
2 


) exp[ —2(a—a:) 


3m—4 m—2 m—2 oo 


2 


Let us rewrite this equation in a more useful 
form for the special case of a dislocation inside 
the plane face of an infinite block of material. 
We find that, to a good approximation, 
mG)? 2R 
= ———— log——+surface terms, (28) 
ro 

where R is the distance of the dislocation from 
the surface of the cylinder. The first term is the 
leading term in this expression. This term repre- 
sents the strain energy in the block due to the 
stresses produced by the dislocation and the 
“image’”’ dislocation which is a distance R outside 
the plane face. The use of an “image’’ dose not 
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Fic. 5. Bipolar coordinate system. 
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Fic. 6. Strain energy 
as a function of dis- 
tance from surface. 


PLANE FACE 


“CYLINDER 
| 


| 


ELASTIC ENERGY (ev per atomic PLANE) 


| 


DISTANCE FROM SURFACE (CM) 


completely satisfy the boundary conditions at 
the plane face and the additional stresses used 
produce the “‘surface terms”’ in the strain energy. 
These “surface terms’’ are almost independent 
of the distance of the dislocation from the sur- 
face except at small distances. Because of this 
the first term increases in importance as the 
dislocation moves away from the surface. If we 
consider a very large copper block which has a 
dislocation 510-7 cm inside a plane face we 
find that the surface terms are 26 percent of the 
image term. If the dislocation is 10-* cm from 
the surface, the surface terms amount to 8 per- 
cent of the image term. The force acting on the 
dislocation is obtained by taking the derivative 
of the strain energy with respect to the distance 
from the surface. The result for the case of a 
dislocation inside a plane surface is 


mGr? 1 
————— —+ surface terms. (29) 


The curves drawn in Figs. 6 and 7 give the strain 
energy and the force acting on a dislocation 
as a function of its distance from the surface in a 
copper block. We have considered two cases, 
one in which we have a dislocation near the plane 
face of a semi-infinite block, another in which we 
have a dislocation in a cylinder of radius 10~* 


cm. In the second case the strain energy reaches 
a maximum when the dislocation is at center of 
the cylinder. The strain energy curve is sym- 
metrical about this maximum. Instead of plotting 
the elastic energy of a slab of unit thickness, we 
have given the elastic energy in a slab whose 
thickness is the distance between successive 
atomic planes which are perpendicular to the 
axis of the dislocation. 

We can now decide whether it is energetically 
easier to form dislocations at surfaces or in 
pairs inside the specimen. In order to make a 
comparison we shall suppose that slip of one 
atomic unit has occurred in both cases along a 
distance R in the slip direction. For a dislocation 
near a surface, R is the distance from the center 
of the dislocation to the surface. In the case of 
two dislocations inside an infinite block, R is the 
distance between the dislocations. By comparing 
the curves of Fig. 6 with that of Fig. 2 one finds 
that, in general, it requires about twice as much 
energy to produce a pair of dislocations as it 
takes to form a single dislocation near the surface. 
This does not necessarily mean that most of the 
dislocations are formed at the surface of the 
specimen. One must also remember that there is, 
in general, a much larger volume in which dis- 
locations can be formed inside the material 
than there is near the surface. 
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FORCE ACTING ON DISLOCATION (DYNES PER ATOMIC PLANE) 


DISTANCE FROM SURFACE 


If dislocation theory is correct, then the forces 
investigated in this section and in the previous 
one will be of great importance for any theory of 
annealing. Taylor has shown that a substance 
containing a large number of dislocations is 
hard. However, if we have a specimen which 
initially contains a large number of dislocations 
per cm?, then the image force acting toward the 
surface will cause dislocations to move out of the 
specimen. If we have positive and negative 
dislocations inside the specimen which are on 
the same slip plane, they can attract and an- 
nihilate one another. Unfortunately, the experi- 
mental data in this particular field are not 
sufficient for a quantitative test of the theory. 
Further data on resoftening temperatures for 
single crystals, on the effect of small amounts of 
soluble impurities on the resoftening tempera- 
tures of single crystals, and studies of the effect 
of the size of the single crystal on resoftening 
would be valuable. Measurements of the in- 
ternal friction as a function of the time for 
material having a decrement of about 10~ on 
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small single crystals of various sizes would also 
be of interest. 


THE ENERGY STORED DURING 
WorK HARDENING 


Experimentally it is found that if a metal is 
severely cold worked a certain amount of energy 
is stored in the metal. This energy can be re- 
leased by heating the material. Taylor and 
Quinney® have made measurements on the 
amount of energy stored in this way. In this 
section we should like to see if we can check 
their results using dislocation theory. 

We shall assume that a work hardened ma- 
terial contains many positive and negative 
dislocations which are arranged in a_ two-di- 
mensional lattice. We shall first calculate the 
total energy of a dislocation pair in the rectangu- 
lar lattice shown in Fig. 8. The elastic energy 
required to form a dislocation pair whose distance 


8G. I. Taylor and H. Quinney, Proc. Roy. Soc. 143, 307 
(1934); 163, 157 (1937). 
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Fic, 8. Rectangular lattice and slip planes. 


of separation is R is given by Eq. (18) as 
mG)? R 


2n(m—1) 


3]. (30) 


The energy of a dislocation pair in a lattice is the 
energy required to form an isolated pair minus 
the binding energy of the pair in the lattice. The 
interaction of two dislocations is given by a 
logarithmic potential such as is shown above. 
The above relation includes not only the inter- 
action between the two dislocations, but also 
the self-energy of the dislocations. 

Let us calculate the interaction energy between 
the dislocation pair located at the origin and the 
pair located at x=0, y=2nb (see Fig. 8). We 
shall assume that the members of each pair are 
separated by a distance 6. The energy which the 
system acquires because of the presence of the 
positive dislocation at 0,—6/2 in the field of 
the pair at 0, 2b is 


The energy which the system acquires because 
of the presence of the negative dislocation near 


the origin in the field of the pair at 0, 2nbd is 


mG)? (~—) 
Thus the interaction energy of the two pair of 
dislocations is 


It is evident that W, vanishes when n ap- 
proaches infinity. W, is therefore an interaction 
energy, not a self-energy. Using this expression 
we can calculate the interaction energy between 
the pair at the origin and all other pairs in the 
plane x=0. We find 


mG? « 4n? 
W.= —— loe( ). (32) 


a(m—1) »=1 4n?—1 


We should like to calculate the interaction 
energy between the pair at the origin and the 
dislocations in the plane for which x=+7a. 
In making this calculation we shall use Eq. (30) 
to obtain the interaction energy of two disloca- 
tions. In using Eq. (30) for any two dislocations 
of the lattice we shall assume that @ remains 
constant although R may change. When the 
lattice of dislocations is considered we may there- 
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Fic. 9. Dislocations in rectangular lattice. 


fore vary the distance of separation of the dis- _ has calculated the potential due to a plane lattice 
locations, but not the angles of the lattice. The of positively charged lines. Using his result we 
resulting interaction potential has exactly the find that the interaction energy between the 
same form as the electrostatic interaction of two positive dislocation at x=0, y= —}3b and the 
long, parallel, uniformly charged rods. Madelung® positive dislocations in the plane at x = +7a is 


cos 
2r(m—1)i=1 b (m—1)b 


where K is a constant. Similarly the interaction energy between the negative dislocations in the 
plane x= +7a and the positive dislocation being considered is: 


mG)? ml(—3b—3b)) mG)*?ra 
cos| {+ 
2r(m—1)i=1 1 b (m— 


The total interaction energy of the positive dislocation with dislocations in the plane x= ra is 
Qe-tlralb 


With this result we find the total interaction energy of the positive dislocation with dislocations in 
all planes except the plane having x=0 is 


(33) 


2mGr2 « etlralb 


2W5=——__ (34) 
r=1 m(m—1) r=1 l 
lis odd 
The sum over / can be carried out and we find that Eq. (34) becomes 
mG)? E tog ( 
x(m—1) 1) log 1 ( 


*E. Madelung, Physik. Zeits. 19, 527 (1918). 
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TABLE I. Data used in calculations. TABLE II. Results of calculations. 
NaCl Cu Al Ni Fe Cu Al Ni Fe 
Slip direction (110) (110) (110) (110) (111) Spacing (a) 34.72 32.38 35.75 29.99 
li (Rectangular lattice) 
Modulus of rigidity : 1.872 4.53 2.665 7.69 8.114 Spacing (a) 29.24 27.26 30.18 25.03 
X (dynes/cm?) (Diagonal lattice) 
Poisson’s ratio 0.2053 0.340 0.343 0.309 0.280 Energy stored (cal./g) 0.518 1.013 0.851 0.963 
(Rectangular lattice) 
ro(A) 7.969 6.15 6.90 5.71 5.48 Remey stored (cal-/g) 0.507 1.130 0.832 0.942 
(A) 3.979 2.552 2.855 2.480 2.478 iagonal lattice 
E stored (cal. 0. j 
Density (g/cm!) 2.165 8.93 2.699 8.80 7.865 =e 


We now note that the interaction energy of the negative dislocation at x =0, y= +40 with dislocations 
in all planes except the plane having x = 0 is equal to the interaction energy of the positive dislocation 
as given in Eq. (35). Thus using Eqs. (30), (32), and (35) the total energy required to form a pair of 
dislocations in a rectangular lattice is found to be: 


mG)? b ‘ 3m ( 4n? 1+e-*relb 


This is not yet complete since we must add to this twice the energy required to distort the material 
inside a cylinder of radius 79 around the center of a dislocation. This cannot be calculated from 
elastic theory. However, H. B. Huntington'® has made a calculation of this energy using atomic con- 
cepts. He found the energy inside a cylinder of radius ro to be about 0.8 ev per atomic plane in the 
case of sodium chloride. We can estimate the value of this energy in other materials in the following 
way. It is assumed that the strain energy inside the cylinder of radius ro is proportional to the strain 
energy outside the cylinder in the case of an isolated dislocation. The constant of proportionality is 
taken to be the same for all materials and is determined by using Huntington’s value for the energy 
inside the cylinder in the case of sodium chloride. For a dislocation at the center of a cylinder of 
material of radius 10~* cm, there is about four times as much strain energy outside ro as there is 
inside. 

It will be noticed that the rectangular lattice of dislocations is not stable. The lattice planes parallel 
to the x=0 plane repel one another and the lattice tends to spread out in the x direction. This is not 
a really serious difficulty for two reasons. First, there exists a small periodic potential of atomic 
origin which probably prevents dislocations from moving along the slip planes unless some external 
stress is present. Second, the actual lattice of dislocations cannot be very stable because resoftening 
takes place in some materials at room temperature. Extensive resoftening can be brought about 
in a single crystal of aluminum by heating it to 600°C. The lattice cannot decrease its energy by 
making 6 smaller because dislocations can only move along slip planes, not perpendicular to them 
Thus 0 is determined by the way in which the dislocations are formed. 

We have also made calculations for a lattice which is more stable than the rectangular lattice. 
The arrangement of dislocations in this lattice is shown in Fig. 9. The calculation of the energy 
required to produce a pair of dislocations in such a lattice is performed in exactly the same manner 
as the calculation for the rectangular lattice. The energy of a pair of dislocations in this case is 


b 1 ( n? ) +2 
s odd 


1° H. B. Huntington, Phys. Rev. 59, 942A (1941). 
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The above equation gives a pair energy which is least when a and 6 are small. A decrease in a is 
presumably prevented by the periodic potential mentioned above. 

Calculations of the energy stored in various materials during severe work hardening have been 
made using Eqs. (36) and (37) together with values for the energy inside 79. The experimental informa- 
tion used in these calculations is collected in Table I and the results are given in Table II. We have 
assumed in these calculations that a equals 6. The calculations are made in two ways. First, one can 
determine that separation distance a which will give agreement with the experimental value of the 
energy stored. It is most convenient here to use as the unit of length the distance between nearest 
neighbors measured in the slip direction. The unit of length will thus be different for each substance 
considered. The results of these calculations are given in the first two rows of Table II. These results 
indicate that, in the most thoroughly work hardened state, dislocations are separated by about the 
same number of atomic spacings in different materials. A second way to make the calculation is to 
assume that dislocations are separated by the same number of atomic distances in different work 
hardened materials. The energy stored during work hardening can then be calculated and checked 
against the experimental value. The results of this calculation are given in the third and fourth rows 
of Table II. The experimental values are given in the fifth row of the table. The dislocations were 
assumed to be separated by 34 atomic distances in the rectangular lattice and 29 atomic distances 
in the diagonal lattice. These results are in remarkably good agreement with experiment except in 
the case of iron. The large error found for iron may be due to the fact that this metal has a differ- 
ent structure and a different set of slip planes. We have assumed that the axes of all line disloca- 
tions of the lattice are parallel. This is probably not the case in iron. 

It might be well to point out that the results obtained with the rectangular dislocation lattice 
do not differ markedly from those obtained with the diagonal lattice. It is probable that neither of the 
lattices considered here exists in a work hardened material. We are, therefore, fortunate in finding 
that the results of the calculation do not depend much on the particular dislocation lattice used. 
The energy stored during work hardening is determined mainly by the elastic constants of the 
material, its density, and the value of the unit slip distance. We do not yet understand why disloca- 
tions are separated by the same number of atoms in different work hardened materials. 

The separation distances of dislocations in severely work hardened materials are about 10-* cm 
according to our calculations. For example, for copper in the case of the rectangular lattice a=0.89 
X 10-* cm; in the case of the diagonal lattice a= 0.74 X 10-* cm. This general result agrees with values 
for the density of dislocations in highly work hardened materials found by using other information." 

The author is much indebted to Professor F. Seitz for suggesting several of these problems and for 
many valuable discussions during the progress of the research. 


"W. F. Brown Jr., Phys. Rev. 59, 528 (1941); F. Seitz and T. A. Read, J. App. Phys. 12, 178 (1941). 
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Letters to the Editor 


Further Confirmation of the Montgomery Theory 
of Counter Discharge 


W. E. RAMSEY AND WAYNE L. LEES 


Bartol Research Foundation of the Franklin Institute, 
Swarthmore, Pennsylvania 


July 14, 1941 


HE Montgomery theory! of counter discharge makes 

it possible to predict the dependence of pulse size 
(maximum negative potential change of the counter wire) 
upon wire capacity and counter length under conditions 
where it is safe to assume that the discharge starts simul- 
taneously at all points in the tube. The picture is particu- 
larly interesting if, as we vary counter length and wire 
capacity, the discharge passes from what has been termed 
a “slow” breakdown to a “fast” breakdown (for non-self- 
quenching gas mixtures). It follows from the discharge 
mechanism that for a given cylinder potential and wire 
capacity there is a length /, such that for all /</, the pulse 
size is constant and equal to the over-voltage on the counter. 
For values of 1>/,, V=la/C, where a@ is the charge per unit 
length in the positive ion sheath necessary to quench the 
discharge at the assigned cylinder potential. For /</, the 
counter is breaking multiply in each discharge and the 
succession of these multiple breaks produces what has been 
called a slow breakdown. For />/, the counter is “over- 
shooting.” Here there is only one break in the discharge and 
the counter has been referred to as having a fast breakdown. 
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Fic. 1. Maximum potential change of counter wire and its variation 
with counter length for a series of wire capacities. C expressed in farads 
in the empirical expression, a =5.2 X10™", 8 =1.2 X10-". Low value 
of A due to end effect in short counter. 

Fic. 2. Counter cylinder composed of nine separately insulated 
segments. The number of segments above (here 225 volts above) starting 
potential defines the effective counter length, other segments being at 
a value below the starting potential. Thus by changing the potential of 
the segments the effective counter length may be changed without 
altering capacity of the wire system. Capacity changes were made by 
adding condensers between wire and ground. Voltage swing measured 
using a weakly coupled electrometer tube and grounding counter wire 
with key after each discharge. Counter filled with 9 cm argon-oxygen 
mixture. 
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If the discharge does not start at all points in the tube 
within a short time interval (10-7 second) after it has 
started at one point, the linear charge density will not be 
constant over the counter length and the picture, although 
unaltered in form, is modified somewhat in detail to allow 
for the finite time required for the discharge to propagate 
itself from one point to another. Simple considerations 
combined with direct observations on the time variation 
of counter wire potential’ indicate that this time of propa- 
gation is only important in modifying the picture for small 
C and large /. 

Figure 1 is a family of curves showing the variation of 
pulse size with / for a series of values of C. Here great care 
was taken in the experimental procedure to obtain low 
values of C, variation of | without altering C, and a true 
value of the maximum negative voltage when no charge is 
allowed to leak from the counter wire. Refer to Fig. 2 and 
its caption for a description of the procedure. 

In Fig. 1, ABC is the characteristic for C=5.0uuf, 
ADE for the wire capacity C equal to 7.5uuf, AFG for 
C=10.0uuf, and so on. Examination of this family reveals 
that the lengths AB, AD, AF, AH, AJ (values of J.) are 
proportional to C. The family may be represented satis- 
factorily by the relation V =Constant =over-voltage for 
values of 1, and V=(l/C)[a—£(l/C)] for values of 
Here a@ is the constant charge per unit length referred to 
above and experimentally obtained for small / and large C. 
Thus to a first approximation propagation time reduces 
the average charge per unit length by an amount propor- 
tional to //C. 

The experimental picture obtained agrees with the 
picture presented by the theory, altered only to allow for 
propagation times of the order of 10~? second. The authors 
have been unable to reconcile the results obtained with any 
picture of the discharge differing materially from that 
presented by this theory. 

1C, G. Montgomery and D. D. Montgomery, Phys. Rev. 57, 1030 


(1940). 
2 W. E. Ramsey, Phys. Rev. 57, 1022 (1940). 


New Reactions in Nickel 
R. L. Doran, Purdue University, Lafayette, Indiana 
AND 
W. J. Henperson, National Research Council, Ottawa, Ontario, Canada 
August 7, 1941 


ONOPINSKI has shown! that the half-lives of many 
radioactive nuclei show certain regularities when 
expressed as functions of their atomic numbers. These 
regularities occur in families of nuclei that differ only in 
the number of alpha-particle units that they contain. If the 
logarithm of the period is plotted against the atomic num- 
ber for a given family, the resulting curves may be in- 
terpolated to make predictions about unknown and un- 
certain members of the family. 

Curves of the families containing radioactive isotopes of 
nickel are drawn in Fig. 1. According to the predictions of 
these curves, Ni®? should have a half-life of about 2 minutes 
and Ni®* a half-life of about 36 hours. A period of 36 hours 
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Fic. 1. Log of halt-life 7 as a function of atomic number Z. 


is known in nickel but it has been assigned? to Ni®’ since 
such an activity was not observed with slow neutron or 
deuteron bombardment of nickel. 

In the present work the 36-hour activity was produced 
by a fast neutron reaction on nickel. In addition, two other 
periods of 2 minutes and 2.6 hours were observed. The 
2.6-hour activity has been definitely assigned* to Ni®. Also, 
alpha-particle bombardments of iron showed the 2-minute 
and the 36-hour activities. The activities are then to be 
assigned to Ni*? and Ni®®, since the stable isotopes of iron 
which could produce radioactive nickel by the (alpha, 7) 
process are Fe and Fe®®. Long bombardments of nickel 
by slow neutrons produced a weak 36-hour activity which 
was followed over 4 half-lives. This period was also pro- 
duced by deuteron bombardments of nickel. In both cases 
the activity was a positron activity of 36 hours half-life, 
and must be due to the formation of Ni®® by Ni®8 (n, gamma) 
and (d, p) reactions. 

The establishment of Ni®® to have a half-life of 36+1 
hours leads to the assignment of the 2 minutes and was 
produced by a Fe* (alpha, m) and a Ni®8 (m, 2”) reaction. 
Chemical separations for nickel were made throughout. 

1 Konopinski and Dickson, Phys. Rev. 58, 949 (1940). 


2? Livingood and Seaborg, Phys. Rev. 53, 765 (1938). 
* Heyn, Physica 4, 1224 (1937). 


The Antisymmetrical Interaction in 
Beta-Decay Theory 
CHARLES L. CRITCHFIELD* AND EUGENE P. WIGNER 


Princeton University, Princeton, New Jersey 
August 2, 1941 


HE process of beta-disintegration can be considered 

as a four-particle process, i.e., a process in which four 
particles are created simultaneously. Thus, e.g., the process 
leading to the emission of an electron can be considered as 
the creation of a proton, a neutrino, and an electron in 
positive energy states, and of a neutron in a negative state. 


If one adopts this point of view, it appears natural to 
look for interactions which are either totally symmetric or 
totally antisymmetric in the four particles. We found that 
a relativistically invariant interaction which does not 
contain any derivatives of wave functions! cannot be 
symmetric in the particles. There is, however, a totally 
antisymmetric relativistically invariant interaction which 
can be used instead of the ‘“‘vector interaction” originally 
proposed.! 

In terms of the customary five invariants (heavy par- 
ticle covariants contracted on light particle contravari- 
ants) the antisymmetric interaction is a sum of scalar, 
axial vector and pseudoscalar invariants. An axial vector 
interaction has already been introduced by Gamow and 
Teller? as an alternative to the polar vector in order to 
allow beta-transitions in which Ai=+1 and 0 (but no 
0—0). The antisymmetric interaction thus permits Ai=0, 
+1, including 0-0 and with stronger transitions with 
Ai=0 than in the pure axial vector interaction. There is no 
change of parity in the allowed transitions. 

Certain qualitative aspects of allowed transitions under 
the antisymmetric interaction can be seen at once. There 
is an additional cause for the transitions Ai=0. In con- 
sequence thereof, the theory gives a larger ratio for the 
lifetimes of the light 4n+2 nuclei to the lifetimes of the 
light odd nuclei than the Gamow-Teller theory gave.* 
This increases the discrepancy in the case of He® but 
decreases it in most other cases. However, the theory does 
not help in understanding the great irregularity in the 
lifetimes of the 4n+2 nuclei. On the other hand, the 
theoretical ratio of the number of transitions to excited 
and ground states of the odd nuclei studied by Grénblom* 
is improved in every case. 

Forbidden transitions in this theory are qualitatively 
the same as under the axial vector interaction. First 
forbidden spectra are of two types: one identical with 
allowed spectra in shape and due to the finite velocity of 
the nuclear particles; the other distinctly different in shape 
and arising from the finite radius of the nucleus. The 
change of nuclear spin is always zero for the first type 
(including 0-0) but in the second type the spin may 
change by +1 or +2. All the first forbidden transitions 
incur a change in parity. 

The nature of the departure of forbidden spectra from 
the form of the allowed spectra can be determined as 
follows. An electron with the full amount of energy avail- 
able in a given process has a shorter wave-length than a 
neutrino with the full amount available to it because of the 
finite rest energy of the electron. Thus the forbidden 
spectrum of any degree, if proportional to a power of 
’nuc/Alight part, Will tend to be more intense at high electron 
energies than at low. The obvious effect on the Kurie 
plot is to make it concave toward the axes at high electron 
energies. None of the carefully measured spectra are of 
this character, and, in fact, in most cases the tendency is to 
the opposite curvature. 

Careful measurements on the spectra of Na™ and P® 
by Lawson® yield Kurie plots which are good straight lines 
at the high energy end. According to the present theory it 
would appear that these transitions are either allowed or 
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first forbidden of the first type; but this indication is not 
quite peculiar to our theory. The above argument concern- 
ing the spectrum suggests that the same conclusion would 
be reached in any choice of invariants except, perhaps, in 
certain special cases. This general result may be inter- 
preted to mean that the long-lived elements are not long- 
lived for the most part because of forbidden spin or parity 
change, or it may mean that the application of beta-decay 
theory must be further refined. 

Second forbidden transitions in this theory show no 
change of parity and A’=+2, +3. The shape of these 
Kurie plots is concave to the axis at high energy for all 
matrix elements. 

We are indebted to Professors Konopinski and Uhlen- 
beck for helpful discussions of this subject at the Washing- 
ton Conference and to the former for access to his review 
article on beta-decay before publication. 

* National Research Fellow. 

1 KE. Fermi, Zeits. f. Physik 88, 161 (1934). 

2G. Gamow and E. Teller, Phys. Rev. 49, 895 (1936). 

‘EE. Wigner, Phys. Rev. 56, 526 (1939). 


4B. Grénblom, Phys. Rev. 56, 508 (1939). 
* J. Lawson, Phys. Rev. 56, 131 (1939). 


The Second Maximum of the Rossi Curve 
Leo BrRovussarD AND ALVIN C. GRAVES 
Department of Physics, The University of Texas, Austin, Texas 
August 11, 1941 

INCE Schmeiser and Bothe! reported that the second 

maximum in the Rossi transition curve at about 17 cm 
of lead is due to small angle showers a number of investiga- 
tions have been made on this maximum. The results have 
been somewhat contradictory as to the height of the 
maximum and angular divergence of showers responsible 
for the maximum. In view of this contradiction in data 
taken mainly with counters, it was thought worth while to 
investigate this problem with a cloud chamber. This 
method would eliminate the uncertainty in angle inherent 
in counter measurements as well as the masking effect of 
counter background. 

The apparatus used consisted of a vertical, counter- 
controlled, cloud chamber and a vacuum tube time delay 
control circuit. By means of a mirror system a camera took 
two pictures 54° apart. All lead was placed above the top 
counter so that non-ionizing particles as well as ionizing 
particles might cause an expansion. 

The results obtained are shown in Fig. 1. More than 200 
showers were recorded for each point near the maximum 
so that the probable error for these points is less than 5 
percent. The Rossi curve for the total showers recorded 
has a definite maximum at about 17 cm of lead. The rise in 
the curve (nearly 30 percent) is well above the probable 
error as computed from the number of showers. The second 
curve is for two-particle showers with angular divergence 
greater than 20°. The similarity of the two curves suggests 
that the maximum is mainly due to this component. So few 
narrow angle showers were obtained that it was impossible 
to construct a similar curve for them or to draw any con- 
clusions as to whether or not they exhibit) different 
behavior in the neighborhood of the maximum. .\ curve of 
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Fic. 1, Total intensity and shower intensity of cosmic rays. 


total intensity on a different scale is included for purposes 
of comparison. It was obtained from the photographs used 
in obtaining the shower curves. 


1K. Schmeiser and W. Bothe, Naturwiss. 25, 669 (1937). 


Production of Mesotrons by Ionizing Radiation 


Witsoxn M. 


Kenyon College Cosmic-Ray Trailer Laboratory, Summit of Mount 
Evans, Idaho Springs, Colorado 


August 4, 1941 


N a series of 3881 random expansions of a large Wilson 

cloud chamber containing five horizontal lead plates 
each 1 em thick one photograph was obtained which can 
be interpreted as the production by a heavy particle of a 
shower containing at least four heavy particles and 
possibly more. (Fig. 1). 

The incoming particle must pass through a half-inch 
of cold rolled steel to enter the chamber. If it were an 
electron it is extremely unlikely that it would be un- 
accompanied by shower particles. The energy of the in- 
coming particle must be quite great to produce such a 
shower, and this would almost insure the presence of shower 
particles if it were an electron. 

The short heavy track appearing directly under the top 
plate is either a mesotron or a proton. The track which 


Pia. Wilson photograph. Each lead phite ist 
em thick. Figures .1 and B were taken at 30° to the lett and right o1 
normal. 
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penetrates the second lead plate and then becomes very 
thick and passes out of the front of the chamber may be a 
proton. The upper part of the track appears. slightly 
brighter than the neighboring rays. This would indicate 
that the particle was a proton. On the other hand the ray 
is coming towards the cameras and is the only one which is 
in the very front of the chamber, and the brightness of the 
upper part of the track may be due to these two conditions. 
If this is true the particle is a mesotron. The great increase 
in ionization in passing through the plate is a little larger 
than that usually observed for a proton. 

A ray on the left-hand side passes through two plates 
without producing secondaries and then disappears in the 
back of the chamber. This particle is either a mesotron or 
a proton. .\ ray to the right passes through three lead 
plates without producing secondaries and is either a meso- 
tron or a proton. It passes out of the back of the chamber. 
The ray bounding the right-hand edge of the shower passes 
through one lead plate without producing secondaries and 
passes out of the back of the chamber. This might be an 
electron but is more likely to be a heavy particle. 

The shower might be interpreted as a nuclear disintegra- 
tion where all the fragments were protons of high energy. 
On the other hand most disintegrations of this type send 
particles equally in all directions. The most likely interpre- 
tation seems to be that the incoming particle is a proton of 
high energy and the shower consists of mesotrons. 

This work is supported by the Rumford Fund for the 
Study of Radiation of the American Academy of Arts and 
Sciences, the Penrose Fund of the .\merican Philosophical 
Society, the Fund for \strophysical Research, and the John 
Simon Guggenheim Memorial Foundation. 


* Guggenheim Fellow. 


Cloud-Chamber Photograph of Slow 
Mesotron Pair 


DonaLp J. HuGues* 
University of Sao Paulo, Sao Paulo, Brazil 
August 10, 1941 


URING July, some 5000 photographs have been taken 
with a cloud chamber in a field of 1165 gauss at an 
altitude of 15,500 feet at San Cristobal mine, Peru. As a 
complete discussion of the results obtained may not appear 
for some time, it seems desirable to report now on one 
picture of some interest. The numerical results given here 
are those reported at the cosmic-ray symposium in Rio de 
Janeiro, .\ugust 3-9, and are admittedly rough, there 
having been no opportunity as yet for precise measurement 
of the film. 

Figure 1 shows several shower electrons. a contamination 
a-particle and a positive-negative pair of heavily ionizing 
particles. Stereoscopic examination shows the pair to be 
in the same plane, which is not the plane of the other tracks, 
and proceeding from a point in or near the front glass plate 
of the chamber. 

The right-hand particle (negative) has an Hp value of 
1.03 10° gauss cm. Its ionization (estimated from the 
original negative) is about 4 to 6 times the electronic 
minimum. Corresponding to this variation, the calculated 
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Fic. 1. Slow mesotron pair. 


mesotron mass will range from 160 to 196 electron masses. 
The ionization density can also be estimated from the fact 
that it is very nearly equal to that of the knock-on electron. 
The known range of the latter (hence its velocity) gives an 
ionization value of 4.5. A value equal to this for the 
mesotron results in a mass of 170. The particle cannot be a 
proton, for with the Hp value observed it would have a 
range of only 1.2 cm in the chamber. 

The energy and angle of ejection of the knock-on 
furnishes an additional independent mass determination, 
As it ends in the chamber its range, and hence its energy, 
can be easily obtained, the latter being 34 kev. The angle 
of ejection is estimated as lying in the range 50°-60°, which 
gives a mesotron mass of 189+24. If the incident particle 
were electronic, the knock-on energy would be 320 kev at 
60° and even more for smaller angles—an impossible value. 

The mass of the left-hand (positive) particle cannot be 
fixed as accurately as the other for there is no knock-on 
and the curvature is not as constant. For the present it 
must suffice to say that the ionization and J/p values are 
very close to those of the negative particle and the mass is 
of the same order of magnitude—about 180 electron masses. 

Though later calculations will probably change to some 
extent the quantitative results, it seems reasonable to 
conclude that we have here observed the creation of a 
pair of slow mesotrons (KE~5 Mev) either in the gas or 
glass wall of the cloud chamber. The production of meso- 
trons at this altitude has also been shown by the occurrence 
of related penetrating, non-multiplying rays in cloud 
chambers. Dr. Wollan has observed several such cases 
during the present experiments. The results with the mag- 
netic field, however, furnish in addition information on 
the mass and energy of the particles. About sixteen heavy 
tracks with measurable curvature have been obtained and 
will be reported on later. 

We wish to thank the Cerro de Pasco Corporation for 
their hospitality and the use of many facilities during these 
observations. The Academy of Science of Brazil and the 
faculty of Sao Paulo University have rendered great assis- 
tance to us during our stay in Brazil for which we are 
deeply grateful. 
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